no  A 0 3 0 9 o 0 

Mesoscale  Air-Sea  Interaction  Group 
Technical  Report 


A MODEL  OF  THE  SPATIAL  STRUCTURE  AND  PRODUCTIVITY 
OF  PHYTOPLANKTON  POPULATIONS  DURING  VARIABLE 
UPWELUNG  OFF  THE  COAST  OF  OREGON 


BY 

J.S.  WROBLFWSKI 


Department  of  Oceanography 
Florida  State  University 
Tallahassee,  Florida  32306 


/ 


SEPTEMBER  1976 


Support  from  Office  of  Naval  Research  t - .»  ■••.wu 
and  National  Science  Foundation 


SCCuMI  T v CLAiS'Fl.'ATiON  or  THIS  PAGC  mtin  Dolo  Enl.r.ilJ 


REPORT  DOCUMENTATION  PAGE  * 


I IICf'OftT  N-KHER 


HEAD  INSTRUCTIONS 
. F1EFORE  COMPLETING  FORM 


I GOVT  ACCESSION  NO.I  1.  RECIPIENT  S CATALOG  NUMBER 


i.  TYPE  OF  REPORT  A PERIOD  COVEREO 


4 TITLE  (on!  f tNUltL-  I i.  TYPE  OF  REPOR 

A Model  of  the  Spatial  Structure  and  Productivity-}  Tech.  Report 
of  Phytoplankton  Populations  faring  Variable  Up-  / 
filing  Off  the  Coast  of  Oregon  ' 


J^AUTHORf., 


( PERFORMING  OHO,  REPORT  NUMBER 


RANT  NUMBER/*) 


iXJ.  S./Uroblewski 


..AM  ALtMCNT.  PROJECT,  f ASK 

AREA  A PORK  UNIT  NUMBERS 


» PERFORMING  ORGANIZATION  name  and  AOORESS 

/^Departments  of  Hcteorribcgy  imd  Oceanograohy  s — s.  > — — i 

Florida  State  University  ‘ ( /&  / NR-f583-231 

Tallahassee,  Fla.  32306  V \ / 


«•  CONTROLLING  OFFICE  NAME  AND  ADDRESS  f I II.  REPORT  DATE 

Office  of  Naval  Research  . — , /fiTZ  n/  ) 1976 

Code  481  (TlX  T 

Arlington,  VA  22217  1 I 132 


14  MONITORING  AGENCY  NAME  A AOORESS//I  lllloomt  tn> n Controlling  Oil  If)  IS.  SECURITY  CLASS.  (»l  |M»  tofort) 


IS,  NUMBER  OF  PACKS 


/3£p» 


Unclassified 


IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS,  DISTRIBUTION  STATEMENT  Cl  Ifilo  Kopotl) 


APPROVED  FOP.  FJBLIC  RELEASE  - DISTRIBUTION  0."  THIS  REPORT  IS  UNLIMITED 


17.  DISTRIBUTION  STATEMENT**!  I A.  obotfal  *nl*r.E  In  Hlotk  30,  II  glllotonl  /ram 


l«.  SUPPLEMENTARY  NOTES 


Department  of  Meteorology 
Technical  Report 


(I-  K £Y  WOHOS  (CmIMu*  ofi  II  nt-'MMry  Hultty  bltcJt  nunlM) 


Numerical  model  Biological  productivity 

Coastal  upwelling  Spatial  structure 

Oregon 

Phytoplankton 


20  ABSTRACT  fConllnuo  on  tovofoo  otgo  II  nototomrf  m4  Conti tf  A f Al*ck  kim*A*0 

During  the  season  favorable  for  coastal  upwelling  off  the  western  boundary  of 
continents,  the  local  circulation  is  strongly  influenced  by  occasional  wind 
"events"  of  several  days'  duration.  Variability  in  the  V7ind  stress  affects 
the  rate  of  upwelling  and  ultimately  the  local  biological  productivity. 

To  investigate  the  relationship  between  wind  events  and  primary  production 
off  the  coast  of  Oregon  in  August  1973,  a time-dependent,  numerical  model  of 


DD  IJAIOJ  1473  EDITION  OF  I NOVAS  I*  OBSOLETE 

S/N  010  2*1114*6601 


SECURITY  CLASSIFICATION  OF  THIS  RACE  (Rim  Doto  tn lotod) 


/fol 


-hiimitv  ciASSincmoK  c?  this  p»gc^)iw  d»i« 
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tating gyre  of  the  two-cell  circulation.  Under  variable  winds  downwelling 
is  not  as  prevalent,  enabling  the  plants  to  utilize  the  upwelled  nutrients. 
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Mr.  Joseph  Wroblewski  is  a biological  oceanographer.  This 
dissertation  is  the  product  of  six  years  work  in  which  we  have  attempted  to 
develop  a modem  understanding  of  the  physical  and  biological  processes 
sufficient  to  model  a ocqplex  oceanic  ecosystem.  These  models  are  capable 
of  answering  specific  hypotheses ; validating  sets  of  field  data  and  gener- 
ating new  research  directions  both  for  the  mathematical  development  and, 
more  inportaitly,  for  the  design  of  field  experiments.  The  ecosystem  model 
described  herein  is  over  simplified  purposely  to  permit  a complete  under- 
standing of  the  system. 

This  work  contributes  a few  new  ideas  in  formulating  the  ecosystem 
equations;  namely  the  treatment  of  multiple  nutrient  compartments  and  the 
simulation  of  zooplankton  egestion.  Great  care  has  been  taken  to  describe 
the  finite  difference  equations  and  the  boundary  conditions  vMch  will  permit 
other  investigators  to  utilize  these  schemes. 

It  has  been  my  great  pleasure  to  train  Joseph  Wroblewski.  He  came 
to  Florida  State  University  as  an  eager  young  biologist  with  a desire  to 
build  mathematical  models  of  marine  ecosystems.  He  didn't  know  any  mathematics 
despite  his  formal  education.  In  this  independent  research,  he  demonstrates  a 
great  deal  of  capability  in  and  understanding  of  biology,  physics  and  mathe- 
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During  the  season  favorable  for  coastal  upwelling  off 


the  western  boundary  of  continents,  the  local  circulation 

,, P- 

is  strongly  influenced  by  occasional  wind  "events*' of  sev- 
eral days’  duration.  Variability  in  the  wind  stress  affects 
the  rate  of  upwelling  and  ultimately  the  local  biological 
productivity. 

To  investigate  the  relationship  between  wind  events 
and  primary  production  off  the  coast  of  Oregon  in  August 
1973,  a time-dependent,  numerical  model  of  the  upwelling 
circulation  was  coupled  to  a complex  model  of  primary  and 
secondary  productivity.  Primary  productivity  is  a function 
of  nutrient  concentration,  light  intensity  and  temperature. ^ 
The  model  dependent  variables  (phytoplankton  nitrogen,  zoo- 
plankton nitrogen,  nitrate,  ammonia  and  detrital  nitrogen) 
are  advected  by  a flow  which  is  influenced  by  bottom  topo- 
graphy and  a variable  wind  stress. 

Advection  by  a two-cell,  upwelling  circulation  is  the 
major  physical  mochanism  leading  to  mesoscale  patchiness  in 
the  phytoplankton,  zooplankton,  detritus  and  nutrient  fields. 
The  numerical  model  predicts  a phytoplankton  and  detritus 
plume  for  which  conriderable  observational  evidence  exists. 
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•Model  predictions  of  daily  primary  production^(78  to  226 
mg  N m"4  dayr  ■L)^during  intermittent  upwclling  arc  paradoxi- 
cally  comparable  to  production  during  strong  upwelling.^ 
When  northerly  winds  arc  strong,  phytoplankton  are "supplied 
with  more  limiting  nutrient  but  experience  a shorter  cuphe- 
tic  zone  residence  time.  The  phytoplankton  are  advcctcd 
offshore  and  down  to  aphotic  zone  depths  by  the  lower,  cy- 
clonically  rotating  gyre  of  the  two-cell  circulation. 

Under  variable  winds  downwclling  is  not  as  prevalent,  en- 
abling  the  plants  to  utilize  the  upwelled  nutrients. 

^ A new  formulation  for  herbivore  egestion  as  a function 
of  food  availability  is  proposed.  Model  herbivore  dynamics 
allow  high  assimilation  efficiency  at  low  grazing  rates. 

In  addition  inhibition  of  nitrate  uptake  by  phytoplankton 
in  the  presence  of  ammonia  is  formulated  in  a manner  con- 
sistent with  available  data.  Parameter  values  of  conven- 
tional, biological  functions  for  growth,  predation  and 
nutrient  regeneration  are  specified  from  either  measurements 
off  Oregon  or  the  literature. 
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in  the  strong  upwelling  case.  The  bottom 
topography  and  the  wind  stress  at  day  4 are 

also  shown.  The  maximum  u and  w velocities 
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in  the  field  are  -2.9  cm  sec  and  1.4x10 
cm  sec"1,  respectively,  (b)  The  gross 
primary  productivity  of  the  water  column 
after  4 days  in  the  strong  upwelling  case, 
(c)  The  spatial  distribution  of  phytoplank- 
ton after  4 days  in  the  strong  upwelling 
case.  Contour  intervals  are  1 ygat  N *-1  , 

19.  Same  as  Fig.  18  but  for  day  7.  The  maxi- 
mum u and  v velocities  in  (a)  are  -5.9  cm 
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Contour  intervals  in  fc)  are  1.4  ygat  N Jl"1 

20.  Same  as  Fig.  18  but  for  day  10.  The  maxi- 
mum u and  v velocities  in  (a)  are  -6.1  cm 
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sec  and  5.2x10  cm  sec  , respectively. 

Contour  intervals  in  (c)  are  1.6  ygat  N A-1 
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Figure 

21.  (a)  Zooplankton  field  after  10  days  in  the 

strong  upwclling  case.  Contour  intervals 
arc  1.8  ygat  N A’* . (b)  Nitrate  field 

after  10  days  in  the  strong  upwclling  case. 
Contour  intervals  arc  2 ygat  NO-j  A~*  . . . 

22.  (a)  Ammonia  field  after  10  days  in  the 

strong  upwclling  case.  Contour  intervals 
arc  0.6  ygat  NH^  A~*.  (b)  Spatial  distri- 

bution ot  detritus  after  10  days  in  tho 
strong  upwclling  case.  Contour  intervals 
arc  0.12  ygat  N A"*  

23.  Same  as  Fig.  18  but  for  day  15.  The  maxi- 
mum u and  w velocities  in  (a)  arc  -3.6  cm 

-1  -2  -1 

sec  and  1.5x10  cm  sec  , respectively. 

Contour  intervals  in  (c)  are  2.0  ygat  N A-1 

24.  Same  as  Fig.  18  but  for  day  20.  The  maxi- 
mum u and  w velocities  in  (a)  arc  -2.8  cm 

-1  -3  -1 

sec  and  7.9x10  cm  sec  , respectively. 

Contour  intervals  in  (c)  are  1.6  ygat  N A-1 

25.  Same  as  Fig.  21  but  for  day  20.  Contour 
intervals  in  (a)  are  1.6  ygat  N A-1  and 

in  (b)  are  2 ygat  NOj  A’1  

26.  Same  as  Fig.  22  but  for  day  20.  Contour 

intervals  in  (a)  are  0.8  ygat  NH^  A"1  and 
in  (b)  are  0.1  ygat  N A 

27.  The  north-south  component  of  the  wind  stress 

near  the  coast  for  the  strong  upwelling  case 
The  data  points  show  the  maximum  daily  gross 
primary  production  calculated  for  each  model 
day 
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1.  INTRODUCTION 


The  mechanisms  controlling  phytoplankton  patchiness 
in  the  ocean  arc  poorly  understood.  By  "patchiness11  the 
author  refers  to  the  three-dimensional  structure  of  the 
phytoplankton  population,  usually  measured  as  concentra- 
tion of  biomass. 

It  has  long  been  suspected  that  the  spatial  hetero- 
geneity of  oceanic  plankton  is  strongly  related  to  the 
variability  of  the  physical  environment  (Bainbridge,  1957; 
Cassie,  1963).  Any  attempt  to  parameterize  the  physical 
environment  in  less  than  its  full  complexity  results  in 
some  loss  of  understanding  of  the  interaction  of  physical 
and  biological  processes  in  creating  phytoplankton  patch- 
iness. Yet,  parameterization  to  some  extent  is  necessary, 
as  one  cannot  look  at  all  relevant  spatial  and  temporal 
scales  simultaneously. 

Twenty  years  elapsed  between  the  first  dynamical 
investigation  of  patchiness  (Kierstead  and  Slobodkin,  1953) 
and  significant,  subsequent  progress  on  the  theory  of  the 
phenomenon  (Platt,  1972).  Theory  on  the  distribution  of 
phytoplankton  in  the  vertical  dimension  has  progressed 
steadily  (for  a review,  see  Patten,  1968),  in  part  due  to 
the  ease  of  sampling  in  the  vertical  relative  to  areal 
sampling.  Only  recent  development  of  instrumentation  and 
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methodology  (Walsh,  1972;  Denman  and  Platt,  197S;  Powell 
ct_  ajU , 1975)  which  enables  quasi-synoptic  sampling  of  the 
phytoplankton  distribution,  has  made  it  possible  to  field 
test  theories  of  patchiness.  Yet,  even  the  most  recent 
theoretical  papers  on  the  spatial  structure  of  phytoplankton 
populations  (Criminale  and  Winter,  1974;  Kamykowski,  1974; 
Platt  and  Denman,  1975;  Dubois,  1975;  Wroblcwski  and 
O’Brien,  1976)  consider  the  physical  advectivc  and  diffu- 
sive processes  to  be  homogeneous  in  the  horizontal  plane. 

The  few  exceptions  (Walsh  and  Dugdale,  1971;  Walsh,  1975; 
O'Brien  and  Wroblcwski,  1976)  have  been  complex  numerical 
models  with  no  analytical  solution  possible. 

biological  processes  regulated  by  the  physical  and 
chemical  environment  can  interact  with  physical  transport 
processes  to  create  large  scale  features  called  plumes  and 
tongues.  The  spatial  structure  of  phytoplankton  on  scales 
below  several  hundred  meters  is  largely  controlled  by 
turbulence  (Platt,  1972;  Denman  and  Platt,  1975;  Powell  e_t 
al. , 1975).  This  small  scale  variability  is  linked  to 
larger  scale  structures  through  the  continuous  turbulent 
dissipation  of  the  latter  (Nihoul,  1975;  Platt,  1975). 
Differences  in  the  physiological  character  of  phytoplankton 
populations  and  the  existence  of  microenvironments  can  also 
lead  to  Small  scale  patchiness  (Platt  and  Filion,  1973; 

Platt  and  Subba  Rao,  1970;  Richerson  el:  al_. , 1970). 

Above  the  1 km  length  scale  phytoplankton  biomass  no 
longer  behaves  merely  as  a passive  scalar  subject  to 
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turbulent  transport.  Advection  becomes  increasingly  im- 
portant relative  to  diffusion  and  biological  processes  in 
determining  the  spatial  configuration  of  the  phytoplankton 
biomass  (O'Brien  and  Wroblcwski,  1973;  Denman  and  Platt, 
1976). 

As  often  observed  length-scale  of  phytoplankton  patch- 
iness in  coastal  upwclling  areas  is  5 to  10  km  (Beers  c_t 
al. , 1971;  Walsh  ot  al^ , 1974  ; Kelley  ct^  al^ , 1975).  Most 
biological  processes  relevant  to  phytoplankton  dynamics 
occur  within  the  euphotic  zone.  The  time-scales  of  interest 
here  range  from  hours  (diel  periodicity  is  quite  evident  in 
phytoplankton  standing  crops)  to  weeks,  i.e.,  several  times 
the  fundamental  time  scale  of  the  mid-latitude  upwelling 
event  (Buyer  and  Pattulo,  1972).  The  formation  and  dissipa- 
tion of  plumes  and  tongues,  i.e.,  water  masses  rich  in 
phytoplankton  biomass  with  predominant  vertical  and  horizon- 
tal extensions  respectively'1',  is  encompassed  within  these 
temporal  and  spatial  scales  (Beers  £t  al^  , 1971;  Walsh  and 
Dugdale,  1971). 

r : 

The  United  Nations'  Scientific  Committee  on  Ocean  Re- 
search (SCOR)  Working  Group  36  during  their  1974  meeting  in 
Kiel,  Germany,  provided  new  guidelines  for  terminology  used 
to  describe  spatial  distributions  of  biological  and  chemical 
variables  in  upwelling  regions.  The  Group  recommended  use 
of  the  term  "tongue"  to  refer  to  spatial  features  with  a pre- 
dominant horizontal  length  scale,  such  as  a shallow  coastal 
bloom  of  phytoplankton  with  a seaward  protuberance  which 
remains  contiguous  with  the  shoreline,  "Plumes"  should  refer 
to  coastal  blooms  with  a significant,  offshore  directed,  ver- 
tical extension.  Thus  the  nutrient  depleted,  thin  layer  of 
low  salinity  water  which  is  discharged  from  the  mouth  of  the 
Columbia  River  and  overlies  the  coastal  water  off  Washington 
and  Oregon  should  be  referred  to  as  the  Columbia  River 
"tongue." 
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During  the  summers  of  1972  and  1973  the  Coastal  Upwel- 
ling  Experiment  (CUE)  was  conducted  off  the  Oregon  coast. 

CUE  was  part  of  the  Coastal  Upwclling  Ecosystems  Analysis 
Program  sponsored  by  IDOE.  Oceanographic  data  collection 
efforts  during  CUE  wore  carried  out  by  four  ships;  the  R/V 
YAQUINA  and  R/V  CAYUSE  from  Oregon  State  University,  the 
National  Oco  inic  and  Atmospheric  Administration  ship  R/V 
OCEANOGRAPHER,  and  the  R/V  T.G.  THOMPSON  of  the  University 
of  Washington. 

Continuous  current  measurements  were  made  by  three 
types  of  buoy  arrays.  An  aircraft  from  the  National  Center 
for  Atmospheric  Research  made  remote  measurements  of  sea 
surface  temperature  and  color,  in  addition  to  collecting 
meteorological  data.  The  biological  component  of  CUE  con- 
sisted of  sampling  surveys  of  the  primary  and  secondary 
production,  coordinated  by  Drs.  L.  F.  Small  and  C.  B.  Miller 
of  Oregon  State  University. 

During  CUE- I in  1972,  it  became  evident  that  the  local 
upwelling  circulation  was  significantly  affected  by  wind 
"events"  of  3 or  more  days'  duration.  Figure  1 is  a pro- 
gressive vector  diagram  (PVD)  of  hourly  wind  data  recorded 
during  the  summer  of  1973  by  an  anemometer  located  on  a 
jetty  off  Newport,  Oregon.  The  plot  is  constructed  by 
placing  vectors  of  the  wind  speed  and  direction  head  to 
tail.  The  PVD  indicates  the  winds  were  predominantly  from 
the  north  during  the  summer  months.  Two  wind  events  in 
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July  ure  shown.  The  first  begins  July  8 and  lasts  7 days, 
and  the  second  begins  July  21  and  endures  an  exceptionally, 
long  20  days. 

Thi  japer  is  an  attempt  to  determine  whether  our  cur- 
rent understanding  of  coastal  upwelling  circulation  and 
marine  biological  processes  can  be  combined  into  a dynamical 
explanation  of  mesoscaic  phytoplankton  patchiness.  Non- 
linear equations  for  the  distribution  of  phytoplankton, 
herbivores,  detritus  and  the  nutrient's,  nitrate  and  ammonia 
in  a transverse  plane  normal  to  the  Oregon  coast  are  solved 
numerically  for  both  strong  and  intermittent  upwelling  con- 
ditions. Daily  primary  production  of  the  water  column  is 
calculated  at  1 km  intervals  within  50  km  of  the  coast 
during  20  days  of  strong  upwelling.  Daily  primary  produc- 
tion is  also  computed  for  the  first  20  days  of  August  1973 
when  variable  winds  induce  intermittent  upwelling. 

The  ecosystem  dynamics  pertinent  to  patchiness  are 
formulated  in  Section  2.  These  dynamics  are  investigated 
first  without  spatial  dependence  and  subsequently  with  ver- 
tical and  horizontal  dependence.  The  steady  state  vertical 
solutions  are  utilized  as  initial  conditions  for  the  time- 
dependent,  two-dimensional  siodel  of  phytoplankton  patchiness 
off  Oregon.  The  numerical  scheme,  boundary  conditions,  and 
physical  dynamics  of  the  latter  model  are  carefully  dis- 
cussed in  Section  3.  Important  questions  concerning  poorly 
known,  biological  parameter  values  and  controversial 
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process  formulations  are  discussed  in  viow  of  an  empirical 
sensitivity  analysis  in  Section  4.  Model  water  column 
productivities  and  dependent  variable  distributions  arc 
compared  to  observations.  Finally,  the  ability  of  the 
model  to  predict  primary  production  and  phytoplankton  patch 
iness  during  coastal  upwelling  off  Oregon  is  evaluated. 


2.  MODEL  FORMULATION  - THE  BIOLOGICAL  DYNAMICS 

In  this  section  the  rationale  for  the  biological  and 
chemical  dynamics  included  in  the  model  is  presented. 

Proper  formulation  of  the  problem  is  the  most  difficult  and 
important  task  in  modeling  research.  The  ability  to  repro- 
duce and  explain  the  behavior  of  oceanic  ecosystems  requires 
the  inclusion  of  the  important  operating  mechanisms  in  the 
model  equations. 

2.1  The  general  equation  for  mesoscale  phytoplankton 
patchiness 

The  general  equation  which  describes  the  distribution 
of  a non- conservative  variable  in  the  sea,  e.g.,  phytoplank- 
ton bioma*  , V , s 

+ V VP  7 • (KVP)  * biological  dynamics  (1) 

where  t is  time,  V represents  the  horizontal  and  vertical 
water  velocities,  and  K is  the  coefficient  of  eddy  diffu- 
sivity.  The  first  term  is  the  local  change,  in  P.  The 
divergence  represents  advection  of  P,  the  third  term  repre- 
sents turbulent  mixing,  and  "biological  dynamics"  refers  to 
the  biological  processes  affecting  the  local  change  in  P. 

Three  fundamental  assumptions  are  necessary.  First, 
the  velocity  field  is  assumed  to  be  nondivergent,  i.e., 

7 • V « 0.  This  is  a requirement  for  conservation  of  mass. 
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Second,  the  horizontal  and  vertical  coefficients  of  eddy 
diffusivity  are  assumed  constant.  Third,  all  derivatives 
in  the  longshore  direction  arc  neglected.  The  region  of 
'the  Oregon  coast  chosen  for  the  major  field  experiment  CUE 
is  an  area  where  this  assumption  is  moro  likely  to  be  valid 
than  other  upwclling  regions  currently  under  study.  If  one 
chooses  a coordinate  system  in  which  y is  in  the  longshore 
direction,  x is  positive  towards  the  coast,  and  z is  posi- 
tive downwards,  (1)  can  be  rewritten, 


iH  + u iP.  + w in 

at  u ax  w az 


- K 


a2P 


K. 


a2p 


Bz 


biological  dynamics 

(2) 


The  horizontal  velocity,  u,  is  assumed  positive  towards  the 
coast,  and  the  vertical  velocity,  w,  is  positive  upwards. 

Platt  and  Denman  (1975)  examined  the  relative  magnitude 
of  the  physical  and  biological  terms  in  an  equation  for  the 
mesoscale  distribution  of  phytoplankton  in  the  ocean,  simi- 
lar to  (2).  They  found  that  given  the  appropriate  condi- 
tions, any  one  of  the  probable  biological  and  physical 
processes  involved  (phytoplankton  growth,  phytoplankton  cell 
sinking,  herbivore  grazing,  advection  or  diffusion)  can 
dominate  the  equation.  Thus  it  is  important  to  specify  as 
exactly  as  possible  the  physical  and  biological  dynamics 
occurring  in  the  region  of  study. 


2.2  Biological  processes  included  in  the  model 

In  the  hierarchy  of  ecosystem  modeling  (Dugdale,  1975) 
the  (x,z,t)  simulation  model  presented  here  could  be 
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classified  as  a productivity  model  with  a detailed  conside- 
ration of  physical  influences  in  both  the  horizontal  and 
vertical  dimensions.  Indeed,  ond  of  the  major  concerns  of 
this  model  is  to  evaluate  the  primary  productivity,  or  the 
rate  at  which  phytoplankton  biomass  increases,  under  differ- 
ent upwelling  conditions.  The  author  attempts  prediction 
of  variances  in  phytoplankton  biomass  as  well  as  simulation 
of  the  standing  crop.  Toward  this  end  one  must  first  iden- 
tify and  properly  formulate  the  biological  processes  influ- 
encing primary  production. 

Neglecting  spatial  terms  for  the  moment,  the  local 
change  in  the  biological  variables  can  be  described  as 


* uptake  of  NOj  and  NH4 
by  growing  P 

- lysis  of  senescent  P cells 
~ * ingestion  of  P by  Z 

0 w 

- metabolic  excretion  by  Z 


- grazing  upon  P by  Z 

C3) 

- egestion  of  fecal 

pellets  (4) 


« fecal  pellet  production 
- bacterial  mineralization 
of  detrital  nitrogen 

3NH4 

3~t  ■ * bacterial  decomposition 
of  detritus 

- bacterial  oxidation  of 

nh4 

3N°3 

—g-jr-  * bacterial  oxidation  of 

nh4 


+ phytoplankton  cell 
lysis  (5) 

+ metabolic  excretion 
by  Z (6) 

- uptake  of  NH4  by  P 


- uptake  of  NOj  by  P 


(7) 
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where  P is  phytoplankton  nitrogen,  Z is  zooplankton  nitro- 
gen, D is  detrital  nitrogen,  NO-  is  dissolved  nitrate  and 

2 

nitrite,  and  NH^  is  dissolved  ammonia  . All  biotic  compo- 
nents arc  expressed  in  units  of  concentration  of  the  limit- 
ing nutrient,  nitrogen.  Biological  rates  arc  expressed  in 
terms  of  nitrogen  turnover  time. 

Park  (1967)  measured  a preformed  nitrate/phosphatc 
atomic  ratio  off  Oregon  of  7:1.  Observing  the  normal  assimi- 
lation ratio  is  16:1,  Park  suggested  nitrogen  as  the  limiting 
nutrient  off  Oregon.  Field  sampling  often  showed  silicate 
and  phosphate  present  in  the  cuphotic  zone  where  nitrate  and 
ammonia  were  depleted  (Ball,  1970;  Atlas,  1973). 

In  most  food  chain  models,  marine  species  with  similar 
feeding  habits  arc  assigned  to  one  trophic  level.  The  indi- 
vidual dynamics  of  species  are  lost  in  this  aggregation, 
except  where  one  organism  dominates  the  trophic  level.  The 
phytoplankton  community  off  the  Oregon  coast  during  the 
upwelling  season  consists  mainly  of  the  diatoms  Skeletonema 
costatum,  Chactoceros  spp.,  Rhizosolenia  spp.  and  Thalassio- 
sira  spp.,  with  dinoflagellates  present  but  less  abundant 
(Anderson,  1972;  Menzies,  personal  doramunication) . The 
zooplankton  over  the  Oregon  continental  shelf  in  summer  are 
mostly  copepoda.  The  copepods  with  the  highest  average  bio- 
mass are  Acartia  clausii , A.  longiremis , Pseudocalanus 
minutus , Caianus  finmarchicus , and  C.  plumchrus  (Peterson, 
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1972;  Myers,  1975).  Over  the  shelf  break  the  pelagic 
species,  Calanus  pacificus  and  Euphausia  pacifica,  become 
important  (Smiles  and  Pcarcy,  1971).  The  collective  physi- 
ology and  behavior  of  these  organisms  specify  the  plankton 
dynamics  included  in  this  model.  For  example,  the  zooplank- 
ton species  over  the  Oregon  shelf  exhibit  little  dicl 
vertical  migration  (Peterson,  1972),  and  thus  this  behavior 
is  not  simulated. 

Consideration  of  ecosystem  dynamics  is  limited  to  two 
trophic  levels  and  a slowly  regenerating  detritus  component 
(Fig.  2).  Carnivore  biomass  or  predation  is  not  considered. 
Reproduction  or  natural  death  of  zooplankton  is  not  dealt 
with  on  the  short  time  scale  of  concern  here  (days  to  weeks), 
although  the  herbivores  can  increase  in  biomass  by  assimila- 
tion of  phytoplankton  nitrogen.  Multiple  nutrient  limita- 
tion of  phytoplankton  growth  is  restricted  to  the  dissolved 
nutrients,  N03  and  NU^.  It  is  assumed  extracellular  excre- 
tion of  nitrogen  by  growing  plants  is  negligible.  Plant 
growth  inhibition  or  enhancement  by  trace  elements  and 
chelation  effects  (Johnston,  1964;  Barber  and  Ryther,  1969) 
are  not  considered. 

2.2a  The  Phytoplankton  equation 

Dugdalc  (1967),  Eppley  and  Coatsworth  (1968),  Maclsaac 
and  Duy^ale  (i960) , and  Caperon  and  Meyer  (1972b)  have  demon- 
strat  that  uptake  rates  of  nitrate  and  ammonia  by  marine 
phytopla.lkton  can  be  expressed  as  hyperbolic  functions  of 


Fig.  2.  Biological  and  physical  processes  included  in  the  spatial  model.  The 
biotic  components  include  phytoplankton,  herbivores,  detritus,  and  the  dissolved  nutrients 
nitrate  and  ammonia.  Carnivore  dynamics  are  ignored  in  the  model,  although  the  figure 
shovs  the  probable  role  of  carnivores  in  the  ecosystem.  Thin  arrows  denote  pathways 
of  nitrogen  between  model  components.  Heavy  arrows  indicate  spatial  fluxes. 
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nutrient  concentration  when  that  nutrient  limits  growth. 
The  Michaclis -Menton  formulation  describing  these  uptake 
kinetics  is 


where  V is  the  uptake  rate  (time'*)  of  nutrient  N (concentra- 
tion), Vm  is  the  maximum  uptake  rate,  and  ku  is  the  Michaclis 
or  half-saturation  constant.  The  concentration  ky  supports 
half  the  maximum  uptake  rate. 

Epplcy  and  Thomas  (1969)  and  Caperon  and  Meyer  (1972a) 
recognized  that  nutriont  uptake  and  cell  growth  are  indi- 
rectly related,  as  uptake  can  bo  separated  in  time  from  cell 
division.  Epplcy  and  Thomas  have  suggested  that  since  cell 
growth,  p,  is  a function  of  cellular  content  of  limiting 
nutrient,  growth  is  also  a hyperbolic  function  of  dissolved 
nutrient  concentration 


Only  if  the  half-saturation  constants  for  uptake,  ku,  and 
growth,  k , are  equal  is  nutrient  uptake  equivalent  to  cell 

o 

growth.  The  author  makes  this  assumption  here,  a not  inde- 
fensible position  when  modeling  phytoplankton  growth  in 
upwelled  waters. 

Phytoplankton  cells  preferentially  take  up  ammonia 
over  nitrate.  Indeed,  the  presence  of  ammonia  inhibits  the 
activity  of  the  enzyme  nitrate  reductase  essential  to  the 
uptake  kinetics  (Packard  and  Blasco,  1974)  and  acts  by 
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reducing  Vm  (N03)  (Dugdale  and  Maclsaac,  1971;  Walsh  and 
Dugdalc,  1972). 

Measurements  of  Vm  (N03)  versus  ammonia  concentration 
in  the  upwclling  regions  of  Peru  and  Northwest  Africa  show 
a wide  scatter  of  data  points  (Dugdale  and  Maclsaac,  1971; 
Maclsaac  et  al.,  1974).  The  data  presented  in  Walsh  and 
Dugdale  (1972)  and  in  Packard  and  Blasco  (1974)  suggest  an 
exponential  rather  than  a linear  decrease  in  Vm  (NO-j) , with 
increasing  ammonia  concentration. 

To  simulate  suppression  of  nitrate  uptake  by  ammonia, 

.UIMU 

Vm  (N03)  is  multiplied  by  the  exponential,  e 4.  The  con- 
centration of  ammonia  where  uptake  of  nitrate  falls  to 
approximately  one  third  Vm  (N03)  is  V"1.  Figure  3 shows  the 
exponential  reduction  in  Vm  (N03)  with  increasing  concentra- 
tion of  ammonia  found  by  Walsh  and  Dugdale  (1972). 

To  be  mechanistically  correct,  perhaps  the  inhibition 
phenomenon  should  be  modelled  as  a competitive  inhibition 
reaction  with  a threshold  effect  (Dugdale,  personal  communi- 
cation). However,  this  formulation  requires  the  specifica- 
tion of  the  dissociation  constant  for  the  nitrate  reductase- 
ammonia  complex  (White  et_  al_. , 1968)  which  is  difficult  to 
measure.  To  avoid  this  complex  formulation,  nitrate  inhibi- 
tion is  modelled  in  an  empirical  rather  than  a mechanistic 
manner. 

In  Fig.  4 the  theoretical  uptake  rates  of  nitrate  and 
ammonia  are  shown  for  increasing  concentrations  of  N03  and 
NH^.  Total  nitrogen  uptake  by  phytoplankton  is  given  by 
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where  it  has  been  justifiably  assumed  Vm  (NO^)  “ Vm  (NH^) 
and  the  half-saturation  constants,  ku,  for  nitrate  and 
ammonia  are  equal  (Eppley  e£  al . , 1969;  Maclsaac  and 
Dugdale,  1969). 

The  loss  of  nitrogen  from  the  phytoplankton  population 
by  cell  autolysis  is  represented  by  a linear  loss  term,  -HP, 
although  the  process  is  a complex  function  of  physiological 
stress.  This  term  is  essential  in  properly  modeling  the 
phytoplankton  dynamics  in  the  aphotic  zone  of  the  water 
column. 

The  grazing  function  is  the  Ivlev  (1945)  equation  as 
modified  by  Parsons  ct^  al_.  , (1967)  . 

R ■ ■» 

* 0 ; P < Pt 

where  R is  the  rate  of  ingestion  (hr~*);  RJn  is  the  maximum 
ingestion  rate;  A (conc"^)  is  the  Ivlev  constant  which  modi- 
fies the  rate  of  change  in  ingestion  with  phytoplankton 
concentration,  P;  and  Pt  is  the  threshold  concentration  of 
phytoplankton  at  which  grazing  begins.  Below  the  controver- 
sial threshold  concentration,  the  zooplankton  starve.  The 
values  for  Rm,  A and  Pt  are  species-specific  (Frost,  1974; 
Mullin  et,  al_. , 1975).  The  grazing  rate  as  a function  of 


i =. 


A(P-P*) 


P > P. 
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phytoplankton  concentration  is  piottoo  for  Calanus  pacificus 
in  Fig.  5. 

Upon  substitution  of  these  formulations,  (3)  becomes 


* NO.  -VNH4  N!l4  * 

s 

-A(P-Pt)’ 

0 * 

p - *. 

1 - e 1 

C8) 


where  Z is  zooplankton  biomass  in  terms  of  nitrogen  concen- 
tration. 


2 . 2b  The  herbivore  equation 

Change  in  zooplankton  biomass  is  taken  as  the  difference 
between  ingestion  and  the  sum  of  egestion  and  metabolic  ex- 
cretion. Ingestion  is  calculated  from  the  Ivlev  equation 
discussed  above.  Egestion  rate  as  a function  of  food  avail- 
ability is  computed  from  the  proposed  expression 

T(P-Pt) 

E 


m 


E+Atc*  v’ 
m 


t'-l] 


where  Em(hr~*)  is  the  maximum  egestion  rate;  A (hr  A)  is  the 
egestion  rate  at  the  grazing  threshold,  Pt;  and  T (cone  *) 
determines  the  increase  in  egestion  rate  with  increasing 
phytoplankton  concentration,  P. 

Conover  (1966)  suggests  egestion  is  a constant  propor- 
tion of  food  ingested  for  Calanus  hyperboreus  feeding  over 
a wide  range  of  diatom  food  concentration.  At  low  phyto- 
plankton concentrations,  however,  egestion  may  no  longer  be 
a linear  function  of  ingestion.  The  sigmoid  shaped  eges* ' an 
rate  curves  shown  in  Fig.  5 were  computed  from  the  above 
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proposed  expression.  Egestion  rate  as  a hyperbolic  curve 
following  Conover's  hypothesis  is  also  shown  for  comparison. 

Steele  (1974)  and  Frost  (1974)  recognized  that  effici- 
ency of  assimilation  defined  as  (ingcstion-egcstion)/ingcs- 
tion  may  be  high  when  food  is  scarce,  efficiency  possibly 
decreasing  as  food  concentration  increases.  In  Fig.  6 the 
assimilation  efficiency  is  calculated  for  a range  of  phyto- 
plankton concentration.  Ingestion  rate  is  given  by  the 
Ivlev  grazing  curve  and  egestion  rate  is  expressed  by  the 
sigmoid  curve.  Above  the  grazing  threshold,  where  the 
assimilation  efficiency  is  zero  by  definition,  the  effici- 
ency rapidly  increases  and  then  decreases  to  a minimum  value. 

Metabolic  excretion  of  nitrogen  varies  with  grazing 
activity,  temperature  and  growth  stage.  Nevertheless,  the 
excretion  process  is  expressed  as  a linear  function  of 
zooplankton  biomass,  -TZ,  where  parameter  r is  assigned  a 
value  in  accord  with  laboratory  measured  values  for  marine 
copepods  under  applicable  environmental  conditions.  Equa* 
tion  (4)  can  then  be  written  as 

T(P-P  ) 

E he  z l 

E 7 * FZ 

m * 

(9) 

2 . 2-c  The  detritus  equation 

Detritus  in  this  model  consists  of  egested  copepod 
fecal  pellets  and  ruptured  phytoplankton  cells.  The  term 
for  bacterial  remineralization  of  detrital  nitrogen  into 
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ammonia  is  -<t*D.  Equation  (5)  in  functional  form  becomes 


3D 

W 


a 


E Ac 
m 


T(P-Pt)Z 


E_  + A{oT^‘,rt)  ~ 1] 
m L 


+ HP  " 4>D 


(10) 


2. 2d  The  ammonia  and  nitrate  equations 

The  nutrient  equations  include  the  source  for  uptake 
by  phytoplankton,  the  sink  for  excretion  of  metabolites  by 
herbivores  and  the  Temineralization  of  detritus.  The  equa- 
tion for  ammonia  is 

3NH,  NH.P 

-W1  ' » + rz  - vm  ' nNH4  'n> 

and  that  for  nitrate  is 

3NO-  N0,P  -YNH. 

ST  ' flN1,4  ' Vn, 

Bacterial  oxidation  of  ammonia  into  nitrite  and  subsequently 
into  nitrate  is  expressed  by  the  term  flh’H . . During  the  CUE 
data  collection  effort,  NO2  and  NO^  were  measured  as  total 
NOj.  Thus,  the  model  dees  not  consider  the  nitrite  inter- 
mediate . 

2 . 3 Scaling  of  the  biological  dynamics 

The  biological  equations  (8)  - (12)  contain  explicitly 
the  parameters  Vm,  ku,  T , 1^,  A,  Pt,  E,  Em,  A,  T,  V,  *,  and 
n,  and  implicitly  the  initial  concentrations  P,  Z,  D,  N03 
and  NH^.  By  scaling  the  equations  the  number  of  parameters 
can  be  reduced.  One  nondimens ional  solution  is  then  equi- 
valent to  solving  several  dimensional  cases.  To  transform 
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back  to  dimonsional  units  one  multiplies  the  nondimcnsional 
solution  by  the  scaling  parameters. 

One  can  examine  all  biological  processes  relative  to 
the  doubling  time  of  the  phytoplankton.  If  time,  t,  is 
scaled  by  V^,  parameter  x * tVm,  where  x is  nondimcnsional 
time.  Also  let  P,  Z,  D,  NO-  and  NH^  be  scaled  by  N’t , the 
total  amount  of  nitrogen  (cone)  in  all  biotic  components 
in  the  upwelling  region. 

Using  these  scaling  relationships,  (8)  - (12)  become 
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and  ■ -y-  ■■  * wNH^ 

_ nw3‘  -*NHJ 

e 

where 

P’  = P/Nt 

V * Z/Nt 

D‘ 

- 3/Nt 

NO^  * N03/Nt 

NHJ  « NH4/Nt 

a 

" VN 

* “VVm 

X * r/Vm 

6 

(17) 
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X - ANt 

« “ 5/Vm 
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3 

P*  « Pt/Nt  . 

The  primes  arc  dropped  for  convenience.  All  quantities  in 
(13)  - (17)  arc  nondimcnsional . P,  Z,  D,  NO,  and  Nll^  arc 
all  fractions;  if  multiplied  by  100  they  represent  the  per- 
cent of  Nt  in  that  biotic  component  at  time  t,  c.g.,  a 
standing  crop.  Time  is  also  nondimcnsional.  For  example, 
if  Vm  is  two  doublings  per  day,  then  t ° 2 is  one  dimension- 
al day.  The  reader  must  bear  in  mind  these  scaling  relation 
ships  when  comparing  the  nondimcnsional  model  results  with 
observational  data. 

2.4  Estimation  of  the  biological  parameter  values 

Equations  (13)  - (17)  represent  a time-dependent,  non- 
spatial,  marine  plankton  model.  One  must  next  determine 
the  proper  parameter  values  for  application  of  these  equa- 
tions to  the  Oregon  upwelling  ecosystem. 

Many  of  the  parameters  are  limited  to  a small  range 
after  scaling  by  Nt  and  Vm.  Since  newly  upwelled  waters 
off  Oregon  have  a maximum  total  nitrogen  concentration  of 
approximately  30  pgat  N (mostly  as  dissolved  nitrate), 

Nt  = 30  pgat  N it’1.  There  exists  a considerable  body  of 
literature  on  both  laboratory  and  in  situ  measurement  of 
the  maximum  growth  rate,  V . Eppley  (1972)  has  summarized 
the  existing  data  on  exponentially  growing  diatoms. 


26 


A maximum  growth  rate  of  2.08  doublings  day’1  is  estimated 
for  neritic  diatoms  under  the  environmental  conditions  of 
the  Oregon  upwelling  season. 

One  ratio  which  arises  from  the  scaling  of  concentra- 
tions in  (8)  - (12)  is  a * ku/Nt>  A typical  value  of  the 
half-saturation  constant  for  neritic  diatoms  in  upwelling 
areas  is  1 pgat  N Jl'1  for  both  nitrate  and  ammonia  (Eppley 
c£  al . , 1969;  Maclsaac  and  Dugdale;  1969).  Upon  scaling 

by  Nt,  the  nondimensional  half-saturation  constant,  a,  is 

-I  - 2 

0.03.  Parameter  a typically  ranges  from  10  to  10  for 

most  any  oceanic  area  (O’Brien  and  Wroblewski,  1973). 

Lower  values  of  a correspond  to  phytoplankton  utilizing 

extremely  small  concentrations  of  the  limiting  nutrient. 

A fit  of  the  function,  V/Vm(NOj)  * e"'1,^H4,  by  least 

squares  to  the  data  found  in  Walsh  and  Dugdale  (1972)  gives 

¥ = 1.462  (ugat  NH^/a,)  . Figure  2 plots  the  exponential 

reduction  in  V (NO-)  using  this  value  for  V.  After  scaling 
111 

by  N^,  nondimensional  tp  * 43.86. 

Phytoplankton  physiological  death  is  a significant 
aphotic  zone  process  (Lehman  et  al^. , 1975).  One  can  approxi 
mate  this  loss  rate  of  phytoplankton  nitrogen  in  terms  of 
the  time  necessary  to  reduce  a light-limited  or  nutrient- 
starved  population  to  approximately  one  third  its  initial 
concentration,  PQ.  Under  no  growth  conditions, 

= e~-T  s e’1  when  x * S~1. 

l0 
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The  tine  scale  E"1  is  called  the  e-F^’ding  rate.  If  it  is 
estimated  the  phytoplankton  standing  crop  is  reduced  to 
PqC-1  in  10  days,  then  nondimcnsional  ? n 0.05.  It  is  for- 
tunate that  sensitivity  analysis  (see  Section  4.5)  shows 
the  model  behavior  is  least  influenced  by  this  paramoter, 
for  exact  determination  of  5 in  nature  is  difficult. 

Parameter  6 is  the  ratio  of  the  maximum  hcrbivoi'c 
grazing  rate,  Rm,  to  the  maximum  phytoplankton  growth  rate, 
V . If  8 is  greater  than  one,  the  herbivores  can  reduce 
the  phytoplankton  standing  crop  through  grazing.  For  con- 
tinuous grazing  Calanus  pacif icus , * 0.01  hr'1  (Parsons 

et  al . , 1967)  and  if  Vm  * 0.04  hr"1,  8 * 0.25. 

Assuming  the  herbivore  population  is  dominated  by 
calanoid  copcpods,-  the  species-specific  Ivlev  constant  is 
in  the  range  0.01  to  0.1  (ugst  N/ft)"1  and  the  grazing 
threshold  is  less  than  or  equal  to  2.5  pgat  N ft"1  (Parsons 
£t  al. , 1967).  Upon  scaling,  nondimensional  X * 0.3  to  3.0 
and  P*  <_  0.08.  A value  for  P*  of  zero  is  used  in  all  the 
following  spatial  and  nonspatial  model  solutions. 

The  maximum  egestion  rate  is  chosen  to  simulate  a 
minimum  assimilation  efficiency  of  65%  (Corner  and  Davies, 
1971)  under  superfluous  grazing  conditions,  i.e.,  Em  * 35% 
of  1^.  Thus,  p = 0.35  8.  There  is  currently  no  published 
laboratory  data  on  the  values  of  egestion  parameters  A and 
T.  Based  on  unpublished  work  by  Hirota  and  deduction,  T 
is  chosed  as  0.15  (pgat  N/ft)'1  or  u = 4.5.  For  near 
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complete  assimilation  of  nitrogen  at  the  grazing  threshold, 
A is  7.2xl0"4  hr"1  or  6 - A/Vm  - 0.01. 

Calanus  finmarchicus  grazing  on  an  algal  diet  in  10°C 
waters  of  the  Clyde  Sea  was  estimated  by  Corner  et  al. , 
(1965)  to  excrete  8-11$  of  its  body  nitrogen  per  day.  A 
value  of  101  of  body  nitrogen  per  day  is  a qualified  guess 
for  the  metabolic  excretion  rate  of  herbivores  grazing  in 
the  neritic  upwclling  region  off  Oregon.  Therefore,  non- 
dimensional  y * 0.10. 

Most  of  the  detritus  in  the  model  originates  from 
copcpod  fecal  pellet  production.  Rcdfield  et  al. , (1963) 
have  suggested  that  much  of  the  soluble  nitrogen  in  newly 
formed  fecal  pellets  dissolves  before  the  particle  sinks 
out  of  the  euphonic  zone.  If  it  is  assumed  that  401  of 
the  ammonia  in  fecal  pellets  leaches  out  in  1 day,  then 
parameter  <J>  * 0.5. 

The  rate  of  oxidation  of  ammonia  into  successively 
nitrite  and  nitrate  in  the  ocean  is  a function  of  tempera- 
ture, pressure  and  bacterial  activity  (Von  Brand  et  al , , 
1940).  From  the  data  provided  by  Von  Brand  et  al. , (1937), 
an  e-folding  time  for  oxidation  of  ammonia  to  nitrite- 
nitrate  appears  to  be  about  25  days.  Thus  w * 0.02. 

2.4a  Steady  state  solution  of  the  nonspatial  model 

The  initial  conditions  chosen  for  P,  Z,  NO^,  NH^  and  D 
indicate  the  behavior  of  (13) -(17).  Figure  7 displays  the 
model’s  response  using  the  above  parameter  values  and 
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the  initial  concentrations  P ■ 0.25*  2 = 0.35,  N03  * NH^  ■ 
0.20,  and  D ■ 0.  Ammonia  and  nitrate  arc  successively 
depleted  by  growing  phytoplankton.  The  availability  of 
phytoplankton  results  in  a gradual  increase  of  zooplankton 
biomass.  Grazing  results  in  fecal  pellet  production  and 
the  appearance  of  detritus. 

The  steady  state  eventually  achieved  is  P * 0.245  or 
7.35  ygat  N Jl’1,  Z » 0.685  or  20.55  ygat  N it"1,  NOj  * 0.001 
or  0.03  ygat  N it  1,  NH^  * 0.014  or  0.42  ygat  N it"1  and  D « 
0.055  or  1.65  ygat  N it"1.  These  concentrations  arc  the  spa- 
tially  averaged,  standing  phytoplankton  and  zooplankton 
stocks,  concentrations  of  detritus  and  dissolved  nutrients 
(ammonia  and  nitrate),  where  no  perturbations  of  the  system 
by  external  forcing  functions  exist.  However,  daily  varia- 
tions of  light  and  wind  stress  forcing  of  the  upper  ocean 
prevent  the  biological  dynamics  from  ever  reaching  a steady 
state,  either  in  time  or  space. 

2.5  Formulation  of  environmental  influences  on 

phytoplankton  growth 

In  addition  to  nutrient  concentration,  the  growth  rate 
of  phytoplankton  is  influenced  by  temperature  and  light 
intensity.  The  author  has  chosen  to  model  phytoplankton 
growth  as  a multiplicative  factor  of  nutrients,  light,  and 
temperature.  Walsh  (1975)  regarded  a single  factor  as 
growth  limiting.  However,  Parson  and  Takahashi  (1973)  and 
Platt  et_  al. , (1975)  suggest  that  primary  production  in  the 
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sea  is  regulated  simultaneously  by  several  environmental 
variables . 


2.5a  The  temperature  function 

Where  photosynthesis  is  light  saturated  and  nutrients 
arc  not  limiting,  the  rate  of  plant  growth  may  be  a direct 
function  of  temperature  (Winter  e£  al. , 1975).  Epplcy 
(1972)  found  under  these  conditions  the  specific  growth 
rate  p (doublings  day’1)  could  be  predicted  from  the  empir 
ically  derived  equation 

p » abcT 

where  a is  0.851  doublings  day’1,  b is  the  constant  1.066, 
c is  1°C"1,  and  T is  temperature  in  °C. 

Most  marine  phytoplankton  experience  a suppression 
of  the  growth  rate  above  some  optimum  temperature.  A for- 
mulation proposed  by  Lassiter  and  Kearns  (1973)  well 
describes  the  "growth  vs.  temperature"  or  p vs.  T curve  of 
many  phytoplankton.  Their  equation  is, 
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in  which  TQpt  is  the  optimal  temperature;  Tm  is  the  tempera- 
ture above  which  the  growth  rate  is  zero;  and  q is  the 
constant  which  modifies  the  rate  of  change  in  p with  T, 

Topt  and  Tm*  FiSure  8 shows  the  p vs.  T curve  for  Skeleto- 
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nema  costatum  using  the  above  equation  with  T 
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Fig.  8.  Phytoplankton  growth  rate  as  a function  •» 

of  temperature.  The  solid  line  is  calculated  from  the  I 

equation  of  Lassiter  and  Kearns  (1973) » while  the  dashed 
line  is  calculated  from  Eppley's  (1972)  general  equation. 
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Tm  ■ 50°C  (Jitts  ct,  {U. , 1964;  Curl  and  McLeod,  1961)  and 
q a 0.2°c’*.  Epplcy's  general  equation  for  photoautotro- 
phic  algae  is  shown  as  the  dashed  line  in  Fig.  8.  As  tem- 
peratures are  usually  less  than  TQpt  in  an  upwclling  area, 
Epplcy's  empirical  relationship  is  adequate  for  modelling 
purposes  here. 

The  highest  temperature  observed  within  SO  km  of  the 
coast  of  Oregon  during  strong  upwclling  in  August  1973  was 
14°C.  At  14°C  the  maximum  doubling  time  of  the  phytoplank- 
ton is  2, OS  doublings  day”1  (Fig.  8),  provided  light  and 
nutrient  arc  not  limiting.  Phytoplankton  growth  relative 
to  this  maximum  is  calculated  for  temperatures  less  than 
14°C  by  normalising  Eppley's  equation  by  pm, 

JL  * (1,066)T  - abcT 

where  a * 0,409,  b a 1.066  and  c * 1°C’1,  The  above  equa- 
tion is  used  to  describe  the  influence  of  temperature  upon 
the  growth  rote  of  phytoplankton  in  the  following  models. 

2.5b  The  light  function 

Light  is  another  important  parameter  which  affects  the 
growth  rate  of  phytoplankton.  The  rate  of  photosynthesis 
appears  linearly  related  to  light  at  low  intensities,  be- 
comes constant  at  optimum  intensities,  and  is  suppressed 
at  high  intensities.  This  response  is  described  by  the 
so-called  "photosynthesi;>  vs.  light  intensity"  or  vs.  I 
curve  (Parsons  and  Takahashi,  1973).  Because  the  maximum 
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rate  of  photosynthesis,  P , depends  on  an  optimum  tempera- 
ture, nutrient  and  light  regime,  plots  of  P^  vs.  1 arc 
normalized  by  P^  (Ycntsch  and  Lee,  1966).  A modification 
of  the  formulation  by  Vollcnwcider  (1965) , 

A.  . 3 

Popt  (i  + Ci/is)i:)x/:i (l  + (oi/is)z)n/‘i  ’ 

is  used  to  describe  the  rate  of  phytoplankton  growth  as  a 
function  of  light  intensity,  I.  is  the  irradiance  for 
which  P^  ■ Popt  [2(l+02)n] _1/2.  pQpt  is  the  maximum  photo- 
synthetic rate,  P , multiplied  by  a function  of  G and  n 
(sec  Fee,  1969).  This  complex  formulation  is  adopted  here 
because  of  its  ability  to  fit  P^  vs.  I curves  exhibiting 
photoinhibition.  Figure  9 shows  this  equation  can  reproduce 

the  response  of  diatom  growth  to  light  intensity  if  I * 

- 2 - 1 

0.07  cal  cm  min  , 0 » 0.175  and  n - 4.3.  The  data  points 
in  Fig.  9 are  from  the  P^  vs.  I curve  for  a mixed  culture  of 
the  diatoms  Skeletonema  costatum,  Nitzchia  closterium > Navi- 
cula  sp.  and  Coscinodiscus  excentricus , as  measured  by 
Ryther  (1956)  in  the  laboratory.  Model  experiments  (not 
shown)  demonstrate  primary  production  would  be  greatly  over- 
estimated if  one  used  a hyperbolic  vs.  I curve,  ignoring 
photoinhibition  effects  at  high  light  intensities. 

2.5c  Phytoplankton  growth  as  a function  of 
light,  temperature"'  and  nutrients 

As  nutrient  uptake  is  coupled  (albeit  indirectly)  to 

photosynthesis  (Goering  et_  al_. , 1964*,  Dugdale  and  Goering, 


RELATIVE  PHOTOSYNTHESIS 


LIGHT  INTENSITY  (cal  cm2  min’1) 


Fig.  9.  Relative  photosynthesis  as  a function  of 
light  intensity.  The  solid  line  is  calculated  from  a modification 
of  the  formulation  "by  Yollenveider  (1965).  The  data  points  are 
from  a laboratory  experiment  by  Ryther  (1956). 
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1967)  one  can  equate  relative  photosynthesis  with  relative 
nutrient  uptake,  V/Vm. 

The  relation  of  light,  temperature  and  nutrient  to  the 
growth  rate  of  phytoplankton  in  the  model  is  formulated  as 
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where  all  terms  arc  nondimensional  and  the  units  of  the 
growth  functions  cancel. 

Light  attenuation  with  depth  in  the  ocean  follows  the 
Beers -Lambert  law 


- 2 -1 

where  I is  the  light  intensity  (cal  cm  min  ) at  depth  z, 
Iq  is  the  light  intensity  immediately  below  the  sea  surface, 
and  k is  the  extinction  coefficient  (m’*)  (Parsons  and 
Takashi,  1973).  IQ  varies  with  time  of  day,  cloud  cover „ 
and  season.  The  die!  periodicity  in  IQ  may  be  expressed  as 

I * I sin 
m 


where  I is  the  light  intensity  at  local  apparent  noon, 
d is  the  day length  fraction  of  a day,  and  t is  time.  When 
the  sine  function  becomes  negative,  Iq  is  set  equal  to  zero. 
Periodicity  is  enforced  by  the  modulo. 

Solar  radiation  measurements  made  from  a surface  buoy 
13  km  off  Sand  Lake,  Oregon,  during  CUE  II  indicate  the 


-s 


37 

day  length,  d,  is  13  hi*s  in  August.  The  shortwave  radiation 
penetrating  the  sea  surface  at  local  noon  on  a cloudless 
day,  Im»  is  about  1.25  cal  cm  min  (Reed  and  Halpern, 
197<1).  The  fit  of  a sine  curve  to  this  insolation  data  is 
shown  in  Pig.  10  where  a daylcngth  of  12  hrs  has  been  as- 
sumed for  simplicity  and  scaled  time  begins  at  dawn. 

In  coastal  waters  incident  radiation,  Iq,  is  reduced 
by  a factor  of  two  within  the  first  few  centimeters  of  the 
water  column  as  ultraviolet  and  infrared  radiation  is 
absorbed  (Parsons  and  Takahashi,  1973).  The  penetration  of 
the  photosynthctically  active  400  to  700  nm  wavelength  is 
then  described  by  Beer's  law.  Combining  Beer's  law  with  a 
sine  function,  the  photosynthctically  active  light  intensity 
at  a given  depth  and  time  of  day  is 

I (z , t)  - 0.5  Im  sin  |*  0 KZ  . 

The  extinction  coefficient,  k is  dependent  on  the  absorption 
of  pure  seawater  and  the  scattering  and  absorption  of  sus- 
pended particles.  Lorenzen  (1972)  expressed  k as 

K ■ KW  + V + KdD 

where  k (n>  ) is  the  extinction  coefficient  of  pure  sea- 

n 

2 

water;  tcp  (m  /mg  N)  describes  the  extinction  by  phytoplank- 
ton,  P;  and  (m  /mg  N)  is  the  extinction  by  colored 
dissolved  substances,  tripton,  and  detritus,  D. 

Pjj  vs.  I curves  in  conjunction  with  Beer's  law  have 
long  been  used  to  estimate  the  primary  production  in  the 
ocean  (for  a review,  see  Patten,  j.968)  . However,  the 
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assumptions  of  a vertically  homogeneous  phytoplankton  dis- 
tribution and  a constant  k,  often  made  in  these  calculations, 
is  invalid.  The  self-shading  phenomenon  of  phytoplankton 
necessitates  numerical  integration  of  production  with  depth 
(Fee,  1969). 

If  the  distribution  of  phytoplankton,  P,  with  depth,  z, 
is  known  at  time,  t,  then 


I(z, t)  - 0.5  I sin 


mod  (t,24)' 


tv 


+ Kp  '6  pWdz] 


where  P(z)  is  the  concentration  (mg  N m’J)  of  phytoplankton 
at  depth  z (m) , Kp  (m  /mg  N)  is  the  extinction  coefficient 
per  unit  concentration  of  phytoplankton,  and  kw  (m-1)  is 
the  extinction  coefficient  of  the  seawater  in  the  absence 
of  any  phytoplankton  (Platt  et  al . , 1975). 

Small  and  Curl  (1968)  integrated  l/zQ  /q  P(z)  dz, 
where  Zq  is  the  depth  of  It  incident  light  penetratration 
to  determine  F,  the  average  concentration  of  phytoplankton 
in  the  euphotic  zone  during  the  1S65  Oregon  upwelling 
seai  <n.  Regressing  k upon  F,  they  found 


k = 0.067  + 0.076  P 

-3 

where  F is  expressed  in  mg  chlorophyll  a m . The  value 
for  k , 0.067  m“\  is  higher  than  expected  for  the  absorp- 

W 

tion  of  light  by  pure  seawater  (0.04  m *)  due  to  nonchloro- 
phyllous,  colored,  dissolved  substances  from  the  Columbia 
River  tongue.  Tripton,  or  nonliving  suspended  matter,  from 
local  river  runoff  was  not  as  important  as  these  colored 
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soluble  substances  (Small  and  Curl,  1968).  If  one  assumes 
a chi  a/N  ratio  of  1/8  for  vigorously  growing  phytoplankton 
with  excess  nitrate  in  the  water  (Antia  ct  al. , 1963),  then 
Kp  off  Oregon  during  upwelling  is  0.095  cm  (pgat  N)  , 

The  chi  a/N  ratio  is  quite  variable  in  nature,  weakening 
the  estimate  of  Kp. 

2. Sd  Steady  state  solution  of  the  (z,t)  model 

The  ability  of  the  above  formulations  to  simulate 
plankton  and  nutrient  dynamics  in  a neritic  water  column 
off  Oregon  is  tested  by  comparing  model  solutions  to  obser- 
vational data.  The  initial  conditions  for  the  solutions 
shown  in  Figs.  11  and  i.3  are  the  steady  state  values  for  P, 
Z,  N03,  and  D from  the  nonspatial  model  (Fig.  7).  The 

(z,t)  model  includes  diffusion  of  the  dependent  variables 

2 - 1 

in  the  vertical,  with  a value  for  Ky  of  1 cm  sec  . A tem- 
perature profile  corresponding  to  that  specified  for  a water 
column  50  km  offshore  (see  Section  3.5)  is  assumed. 

Figure  11  displays  the  phytoplankton  vertical  profiles 
which  result  when  plant  growth  varies  with  depth.  The  pro- 
file fluctuates  about  a steady  state  as  incident  radiation 
follows  the  sinusoid  curve  (Fig.  10)  and  photosynthetically 
active  light  penetrating  the  surface  is  attenuated  with 
depth.  A phytoplankton  maximum  occurs  at  a depth  of  21  m. 
Sampling  profiles  taken  50  km  offshore  during  CUE  often 
showed  a chlorophyll  maximum  between  10  and  25  m. 
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NORMALIZED  FRACTION  OF  Nt 


SO. 0-* 


Fig.  XI.  Steady  state  oscillation  in 
the  phytoplankton  profile  of  the  (z,t)  nodel. 
Plant  biomass  is  expressed  as  a nondioensional 
fraction  of  the  total  amount  of  nitrogen  in  the 
upvelling  ecosystem, 
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Fig.  13.  Steady  state  vertical  profiles 
of  the  (z,t)  model's  dependent  variables  (herbi- 
vores, detritus,  nitrate  and  ammonia)  at  midday. 
The  dependent  variables  are  expressed  as  a non- 
diroensional  fraction  of  the  total  amount  of 
nitrogen  in  the  upvelling  ecosystem,  N. . 
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Figure  12  shows  the  light  intensity-depth  profile 
at  midday.  If  light  absorption  by  phytoplankton  wore  not 
considered,  i.c.,  if  ■ 0,  the  1%  Iq  light  intensity 
would  reach  68  m instead  of  33  m.  Suppression  of  phyto- 
plankton rate  by  light  intensities  greater  than  0.13  cal 
-2  -1 

cm  mm  occurs  at  midday  in  the  upper  12  m.  Suppression 
ox  phytoplankton  growth  at  local  noon  in  the  summer  has 
been  observed  to  depths  of  at  least  10  m (Small,  personal 
communication).  Below  IS  m,  the  availability  of  light 
becomes  growth  limiting. 

The  nutrients,  detritus,  and  zooplankton  profiles  at 
midday  arc  shown  in  Fig.  13.  Oscillations  similar  to  that 
shown  in  the  phytoplankton  profile  (Fig.  11)  are  found  in 
the  nitrate  and  ammonia  profiles  but  are  not  depicted.  The 
ammonia  and  nitrate  profiles  indicate  nutrient  depletion  in 
the  upper  20  in  of  the  water  column,  both  nutrients  increas- 
ing in  concentration  below  this  depth  as  uptake  by  phyto- 
plankton is  reduced.  The  zooplankton  maximum  occurs  at  16 
m.  The  zooplankton  standing  stock  quickly  decreases  below 
the  euphotic  zone.  Zooplankton  biomass  is  directly  related 
to  the  concentration  of  its  prey  through  the  Ivlev  grazing 
term.  Detritus,  which  is  assumed  to  sink  at  a rate  of  8 m 
day  1 (see  Section  3.2)  shows  a maximum  at  26  m.  If  a sink 
ing  rate  of  zero  were  assumed,  the  shape  of  the  detritus 
profile  would  more  closely  correspond  to  the  zoof ' xnkton 
profile,  since  egestion  of  herbivore  fecal  pelJ  § is  the 
main  source  of  detritus. 
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2.6  Calculation  of  the  Daily  Gross 

Primary  Production 

Given  tho  phytoplankton  vertical  distribution,  one  can 
calculate  the  rate,  V,  at  which  inorganic  nitrogen  is  incor 
pornted  into  phytoplankton  nitrogen,  P.  Integration  of  tho 
growth  of  phytoplankton  over  depth,  /VP,  where  V is  a func 

M 

tion  of  light,  temperature  and  nutrients  gives  the  gross 
primary  productivity  (mg  h'  m"2  hr"1)  of  the  water  column. 
Models  where  the  functional  expression  for  V is  not  too 
complicated  can  be  integrated  analytically,  provided  the 
vertical  distribution  of  P is  homogeneous,  or  at  best  a 
smooth  function  (Fee,  1970;  Platt  et  al . , 1975).  In  the 
present  model  one  must  resort  to  numerical  integration. 

To  find  the  daily  gross  primary  production  of  the 
water  column  (mg  N m ) , one  integrates  over  depth  and  time 
the  expression, 


VP 

m 


-TNH. 

N°-  e * NIl4  ’ 

+ * *U  + Nll4 


• W, 

[l  ♦ Ci/rsJ2]x/ii [l  ♦ C0l/ls)2ln/i£ 
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where 


1(2, t)  - 0.5  Xm  sin  ^i24? 


-[V+Kp;0  pWdz] 


The  daily  gross  primary  production  calculated  for  the 
phytoplankton  profile  shown  in  Fig.  11  is  100  mg  N m~2. 
Anderson  (1964)  observed  a daily  net  primary  production  of 
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between  0.3  and  1.2  g C m’2  day’1  in  oceanic  waters  50  km 
off  Oregon  during  the  1962  upwelling  soason.  Assuming  a 
C/N  ratio  of  7 (Small  and  Hamburg,  1971),  this  production 
in  terms  of  nitrogen  is  43  to  171  mg  N m day  . It  ap- 
pears the  biological  and  chemical  dynamics  formulated  in 
this  model  can  correctly  simulate  primary  production  off 
Oregon. 


3.  MODEL  FORMULATION  - THE  PHYSICAL  DYNAMICS 


In  this  section  the  physical  dynamics  of  upwolling 
are  formulated  for  incorporation  into  the  spatial,  primary 
productivity  model  in  an  attempt  to  simulate  mososcalc 
phytoplankton  patchiness  observed  during  upwolling  off 
Oregon. 

3 . 1 Wind  forcing  and  bottom  topography 

Thompson  (1974)  has  developed  an  (x,z,t)  model  of  the 
transverse  circulation  off  the  Oregon  coast  for  the  same 
period  in  1973  during  which  on  intensive  biological  sampling 
program  was  conducted.  Thompson's  model  was  run  with  an 
observed  wind  stress  profile  (Fig.  27)  to  simulate  upwolling 
events  which  occurred  in  August  1973.  The  forecast  time- 
dependent  velocity  field  was  used  to  advect  the  biotic  vari- 
ables in  what  will  henceforth  be  referred  to  as  the  "inter- 
mittent upwolling  case."  A "strong  upwelling  case"  was  also 
run  with  a steady  wind  stress  driving  the  upwolling  model. 

Thompson's  circulation  model  delineates  the  position  of 
the  pycnocline,  the  .localities  of  upwolling,  convergences 
and  divergences  all  of  which  are  influenced  by  bottom  topo- 
graphy and  a variable  wind  stress.  Thompson  assumed  a free 
surface  and  the  bottom  topography  shown  in  Fig.  14.  This 
slope  is  a linearized  version  of  the  actual  bottom  topography 
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Fig.  lh.  Bottom  topography  assumed  in  the  (x,z,t) 
upwelling  circulation  model  of  Thompson  (197*0.  The  rectangular 
hatched  region  at  the  upper  right  delineates  the  region  of  the 
(x,z,t)  model  of  phytoplankton  production. 
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off  Oregon.  Bottom  depths  from  tho  coast  to  150  km  offshore 
between  4 4 ° 5 S * N and  44°15'N  were  digitized  from  bathymetry 
charts,  averaged  in  the  longshoro  direction,  and  finally 
linearized  to  obtain  a bottom  slope  characteristic  of  the 
CUK  region.  A bottom  depth  of  50  m at  the  coast  is  assumed 
to  simplify  computations. 

A basic  assumption  of  both  the  physical  and  biological 
spatial  models  is  no  longshore  variation  in  the  coastline, 
bottom  topography  or  velocity  field.  This  assumption  makes 
for  computational  economy.  It  also  simplifies  the  spatial 
dynamics  to  a degree  that  fundamental  features  in  tho  physi- 
cal circulation  and  in  the  biological  and  chemical  fields 
can  be  examined  without  the  confusion  imposed  by  longshore 
variability.  The  importance  of  longshore  variation  in  coast- 
line and  bottom  topography  upon  the  upweliing  circulation 
has  been  discussed  by  Hurlburt  (1974),  Shaeffer  (1974),  and 
Peffley  and  O'Brien  (1976).  However,  the  author  leaves 
consideration  of  longshore  derivatives  in  (2)  to  future, 
more  ambitious,  undertaking. 

All  biological  simulations  are  confined  to  the  upper  50 
ni  of  the  water  column  in  a region  within  50  km  of  the  coast. 
The  modeled  cross  section  is  shown  as  the  thin  hatched  area 
in  Fig.  14.  This  area  is  divided  into  a grid  of  50  by  20 
rectangles  in  the  x and  y directions,  respectively.  The 
dimensions  of  each  grid  box  are  2.5  m in  depth  and  1 km  in 
width.  The  first  biological  grid  point  is  1.25  m below  the 
liea  surface. 
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3.2  Scaling  of  the  physical  dynamics 

The  equation  governing  phytoplankton  patchiness  in  a 
transverse  plane  normal  to  the  Oregon  coast  is  again 

H * u 8 ♦ « I • Kh  0 • Kv  * biological  tor.a  (2) 

3X  d Z 


where  "biological  terms"  refer  to  the  phytoplankton  dynamics 
developed  in  the  previous  section. 

If  one  scales  (2)  so  all  terms  are  nondimcnsicnal,  one 
can  compare  the  relative  influence  of  the  physical  and  bio- 
logical processes  in  determining  the  spatial  distribution 
of  phytoplankton.  Let. 

x * Lx1  u « Uu' 

z c Hz'  w * Ww1 

where  L and  H arc  characteristic  horizontal  and  vertical 
length  scales,  U is  a typical  value  of  the  organized  hori- 
zontal flow,  and  K is  a typical  value  of  the  vertical 
velocity.  Using  the  scaling  relations  put  forth  in  the  pre- 
vious section,  (2)  becomes  nondimensicmal , 


32P* 


scaled  biological  dynamics. 


(18) 


Once  again  the  primes  are  dropped  for  convenience.  The 
scaling  of  the  equations  for  herbivores,  detritus,  ammonia 
and  nitrate  is  similar. 
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Scaling  parameters  L and  !t  can  be  regarded  as  the 
length  scales  within  which  the  pa  'hiness  phenomenon  occurs. 
The  resolution  of  patchiness  is  determined  by  the  choice  of 
Ax  and  Az.  Since  the  author  is  interested  in  resolving  fea- 
tures less  than  50  km  horizontal  length  and  less  than  50  m 
vertical  length,  Ax  and  Az  must  be  several  factors  smaller 

than  these  scales.  Based  on  these  horizontal  and  vortical 

5 2 - 1 

length  scales,  a value  for  of  5x10  cm  sec  and  for  Ky 
2 -1 

of  1 cm  sec  is  assumed  (Okubo,  1971).  The  dimensional 
value  of  all  parameters  used  in  the  (x,z,t)  model  for  both 
the  strong  upwclling  case  and  the  intermittent  upwelling 
case  is  presented  in  Table  1. 

Platt  and  Denman  (1975)  performed  a similar  scaling  of 
(2),  although  they  elected  not  to  divide  through  by  Vm  to 
retain  units  of  frequency.  The  characters tic  time  scale  of 
each  term  is  then  the  reciprocal  of  its  coefficient.  The 
nondimensional  coefficients  in  (18)  define  the  importance 

l 

of  the  adveetive  terms  relative  to  the  diffusive  and  growth 
terms  in  the  equation.  When  U/(LVm)  or  W/(HVm)  >>  1,  advec- 
tion  plays  a dominant  role  in  determining  the  spatial  con- 
figuration of  the  phytoplankton  biomass  (O'Brien  and 
Wroblewski,  1973).  If  one  evaluates  the  magnitude  of  these 
coefficients  using  the  values  presented  in  Table  1,  it 
becomes  evident  that  vertical  advection  is  twice  as  important 
as  horizontal  transport  in  determining  the  spatial  configura- 
tion of  the  phytoplankton.  The  diffusion  terms  are  two 
orders  of  magnitude  smaller  than  the  adveetive  terms. 
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Tabic  1 

Parameter  values  of  the  (x,z,t)  mcsoscale  phytoplankton 

patchiness  model 


Para- 

meter 

Value 

Para- 

meter 

Value 

a 

4.09X10"1 

Vm 

2.4xl0’5  sec’1 

b 

1.07 

ws 

0.01  cm  sec’1 

c 

1.0'C"1 

W 

-2  -1 
2.0x10  cm  sec 

d 

4.32xl04  sec 

Ax 

1.0x10°  cm 

Em 

lxl0“6  sec*1 

Az 

2.5xl02  cm 

H 

5.0x10°  cm 

r 

2.4xl0’6  sec’1 

h 

1.2xio” 5 cal  enf2 

sec’1  A 

2.1xl0"7  sec’1 

2.1xl0”2  cal  cm”2 

sec"1  n 

4.30 

h 

5.0x10°  cm2  sec’1 

8 

0.175 

Kv 

2 - 1 

1.0  cm  sec 

K 

n 

9.5xl0"2  cm2  (ygN)"1 

ku 

1.0  ygN  liter"1 

r 

Kw 

6.7xl0"4  cm’1 

L 

5.0xl06  cm 

A 

0.06  (ygN/liter)’1 

Nt 

30.0  ygN  liter’1 

M 

M 

1.2xl0‘6  sec’1 

pt 

0.0  ygN  liter’1 

T 

0.15  (ygN/liter)  1 

R_ 

m 

2.8x10’°  sec"1 

* 

1. 2x10°  sec’1 

At 

8.6xl02  sec 

y 

1.46  (ygN/liter)"1 

U 

10. 0 cm  sec’1 

n 

4.7xl0'7  sec’1 

S3 


Detritus  (i.c.,  ruptured  phytoplankton  colls  and 
egested  copepod  fecal  pellets  of  floccular  composition)  is 
assumed  to  have  a constant  sinking  rate,  w . The  nutrionts, 

O 

NK^  and  NO-,  are  totally  passive,  i.c.,  wg  * 0.  Phyto- 
plankton cells  may  bo  neutrally  buoyant  in  nutrient  rich 
waters  (Smayda,  1970)  permitting  ono  to  omit  a sinking  term 
in  (2).  The  detritus  sinking  rate,  w , is  sealed  in  the 
same  manner  as  the  vertical  velocity,  w, 


w„  « W w 1 . 
s s 

The  total  derivative  for  detritus  may  then  be  expressed  as, 
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3D 

3X 


+ S. 


cw  + ws>  H 


E a!n  . E A 

11  i? 


m 


sealed  biological  dynamics  (19) 

where  Sr  - U/(LVm),  S2  - W./(HVm),  Eh  - Kh/(L2Vm),  and 


Ev  - V • 

3.3  The  finite  difference  scheme 

Tho  final  nondimens ional  equations  arc  expressed  in 
finite  difference  form  and  solved  for  each  spatial  grid 
point  using  the  scaled  u and  w velocities  derived  from  the 
physical  upwelling  circulation  model.  If  the  vortical 
grid  lines  in  the  x-direction  arc  indexed  by  the  letter,  j; 
the  horizontal  lines  in  the  z-direction,  by  the  letter,  k; 
and  nondimensional  time  is  denoted  by  the  index,  n;  then, 

P?  v * P(ndr,  j£x,  kdz). 

Refer  to  the  stencil  shown  in  Fig.  15.  The  values  of 
the  scaled  u and  w velocities  are  determined  for  the  center 


1 

T 

Fig.  15.  Stencil  showing  grid  point  locations  li 

of  the  u (squares)  and  w (crosses)  velocities  and  the 
dependent  variable  P (circles)  used  in  integrating  T 

the  (x,z ,t)  biological-chemical-physical  model.  The  Ji 

vertical  grid  lines  in  the  x-direction  are  indexed 
by  the  letter,  j,  and  the  horizontal  lines  in  the 

z-dlreccion  by  the  letter,  k.  I 
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of  the  grid  walls  to  increase  computational  stability.  The 
computed  concentration  of  P,  2,  D,  NO-  or  NH^  represents  an 
average  over  the  grid  box.  The  change  in  concentration  of 
a biotic  component  is  dependent  on  the  amount  present  within 
the  grid  box,  the  concentration  in  the  nearest  four  adjacent 
boxes,  the  transport  into  and  out  of  the  box,  and  the  bio- 
logical dynamics  occurring  within  the  box. 

The  finite  differencing  of  the  equations  for  the  bio- 
logical variables  P,  Z,  D,  NO-  and  NH^  incorporates  a leap- 
frog scheme  in  time,  a quadratic  conservative  advective 
scheme  suggested  by  Piacsck  and  Williams  (1970)  and  an 
explicit  scheme  for  the  diffusive  terms.  The  diffusive  and 
biological  terms  arc  lagged  in  time.  For  example,  (18)  is 
rewritten, 


»n+l 

j.k 


pn-l 

3,k 


2 MU 

mN~ 

m 


.*/  -I*1 
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u 


n 

j+l.k 


(P-. 

Il±jt 


i,k+_zLii 


2 At  W 

Az  HV_ 
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where 


rtP?  v+1  ♦ p?  k)  n CP-  k + P?  k-l>  n 
[ }'Y  \ 3,k'  w",kJ 


2&T  Kh  r„n-l  + pn-l 

(ix)2L2v  <-  i+1'k 

2iT  rpn-l 

!jPT  Li'* 


2pi;il 


(Az) 


+ Pn_1 

k+1  Pj  >k-l 


,n-ll 

■j.k] 
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M .k^s 


I Zn-1  - £ 
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When  advcction  and  diffusion  are  neglected,  the  finite  dif- 
ferencing scheme  reduces  to  the  common  Euler  method. 
Conditions  necessary  for  stability  are: 


At  < 


ft  u„ 


where  u,  * Max.  . nr  u 
m j,k,n  LV 


4t  < A* 


At  < 


(4X)  2l\ 


w 

where  wm  B Max.  . _ mr-  w 
m J,k,n  I IV 


At  < 


(42)2I12V„ 


In  this  model  the  advective  criteria  are  the  more  stringent. 
The  stability  criterion  for  the  Euler  method  is  At  < a1 
where  a is  the  largest  biological  rate  parameter.  The  accu- 
racy of  the  finite  differencing  scheme  in  approximating  its 
analogous  continuous  derivative  increases  with  smaller 
values  of  At.  The  time  step,  At,  used  was  0.02  which  con- 
verts to  0.01  days  in  real  time. 
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3 . 4 Boundary  conditions 

No  advoctivo  mass  flux  is  allowed  across  the  solid 
coast  or  the  air-sea  interphase.  At  the  water  column 
bottom  boundary,  the  direction  of  the  Clow  determines  the 
advcctivc  boundary  condition.  IC  water  is  upwolling,  tho 
boundary  values  for  P,  2,  NO-,  NH(j  and  D arc  specified  from 
observational  data.  If  water  is  downwclling,  the  concentra- 
tion of  the  variable  just  inside  the  boundary  determines  the 
value  at  the  boundary.  A one-sided  derivative  is  used  for 
the  advcctivc  terms.  The  offshore  boundary  is  treated  simi- 
larly. No  diffusive  mass  flux  is  allowed  through  any 
boundary. 

These  boundary  conditions  are  such  that  limiting  nutri- 
ent can  become  stored  within  the  model  region  as  either  P, 

Z,  NO-,  NH^  or  D.  Yet,  the  total  amount  of  nutrient  is  con- 
served in  a balance  between  v/hat  is  advocted  into  and  out  of 
the  region  and  what  is  stored  as  standing  stock  or  dissolved 
nutrient. 

Sampling  profiles  of  phytoplankton  nitrogen,  nitrate 
and  ammonia  taken  50  km  seaward  of  the  Oregon  coast  in  1972 
and  1973  closely  approximate  the  steady  state  profiles  of 
these  variables  shown  in  Figs.  11  and  13.  Continuous  pro- 
files of  zooplankton  biomass  and  detrital  nitrogen  at  50  km 
offshore  were  not  taken,  but  the  zooplankton  net  tows  and 
the  particulate  nitrogen  data  do  indicate  the  model's 
steady  state  profiles  for  zooplankton  and  detritus  are 
reasonable.  The  steady  state  values  of  these  variables  at 
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SO  m arc  taken  as  the  boundary  condition  concentrations  in 
upwclling  water  entering  the  model  region  from  below. 

(Table  2). 

3.5  Initial  conditions 

Presently,  there  docs  not  exist  adequate  field  data 
to  specify  the  initial  conditions  of  the  biological  and 
chemical  dependent  variables  for  all  x and  z at  the  onset 
of  the  model  case  runs.  For  lack  of  a better  alternative, 
the  initial  conditions  arc  taken  as  the  steady  state  solu- 
tion of  the  (x,z,t)  primary  productivity  model  in  the 
absence  of  advcction. 

The  temperature  field  is  specified  from  observations 
made  during  August  1973.  Sea  surface  temperatures  were 
found  to  decrease  from  14°C  at  an  offshore  distance  of  50 
km  to  9°C  at  the  coast  during  strong  upwelling.  Tempera- 
tures below  the  surface  mixed  layer  decreased  rapidly  with 
depth,  especially  where  cold  slope  water  moved  onto  the 
continental  shelf.  A polynomial  function  was  fit  to  the 
observed  temperature  data  to  giye  a time-invariant,  smoothed 
temperature  field  evaluated  at  each  grid  location  in  the 
model  region  (Fig.  16a). 

The  influence  of  this  temperature  distribution  on  the 
steady  state,  daily  primary  production  of  the  water  column 
is  shown  in  Fig.  16b.  Fig.  16c  shows  the  initial  conditions 
for  the  phytoplankton  field.  The  phytoplankton  spatial  dis- 
tribution differs  from  the  steady  state  profiles  shown  in 
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Fig.  l6.  (a)  The  spatial  distribution  of  seawater 

temperature  assumed  in  the  (x,z,t)  model  of  phytoplankton 
production.  Contours  are  from  7°C  at  50  m depth  at  the 
coast  to  lJ<°C  at  the  surface  50  km  offshore.  The  contour 
interval  is  0.l4°C.  (b)  The  steady  state  gross  primary 

productivity  of  the  water  column  for  the  (x,z,t)  phyto- 
plankton production  model,  (c)  The  steady  state  phyto- 
plankton spatial  distribution  for  the  (x,z,t)  model  in 
the  absence  of  advection.  Contour  intervals  are  0.8 
ygat  N fl. 
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Fig.  11  only  in  that  plant  growth  as  a function  of  tempera-* 
turo  varies  in  x as  well  as  z (Fig.  16b).  If  the  tempera- 
ture field  (Fig.  16a)  showed  no  variation  with  distance 
offshore,  the  initial  phytoplankton  field  would  appear 
homogeneous  in  x,  although  stratified  in  z. 

At  the  beginning  of  the  modclod  20-day  period  of  the 
strong  upwolling  case,  a spatially  variable  wind  stress  is 
specified  which  linearly  increases  from  zero  at  time  zero 
to  -0.5  dyne  cm  from  the  north  (Fig.  17a  and  b).  The  wind 
stress  remains  at  this  magnitude  from  day  1 to  day  10,  then 
linearly  decreases  to  zero  during  day  11.  It  remains  zero 
for  the  rest  of  the  20-day  period  (Fig.  17b).  The  wind 
stress  used  in  the  physical  circulation  model  for  the  inter- 
mittent upwclling  case  is  derived  from  actual  wind  measure- 
ments recorded  by  on  anemometer  located  on  the  south  jetty 
off  Newport,  Oregon  (Figs.  27  and  28). 


4 . MODEL  RESULTS 


In  this  section  simulations  of  mesoscale  phytoplankton 
patchiness  during  strong  upwclling  and  during  intermittent 
upwelling  in  August  1973  are  presented.  Result  from  both 
model  runs  ore  subsequently  compared  with  observations  made 
during  CUE  and  with  the  literature.  The  goal  is  to  investi- 
gate the  biological,  chemical  and  physical  dynamics  involved, 
rather  than  to  merely  reproduce  the  observed  features  numeri- 
cally. A sensitivity  analysis  of  the  (x,z,t)  model  is  per- 
formed to  ascertain  which  of  the  mechanisms  involved  are  the 
most  important  in  determining  the  spatial  solutions. 

4 . 1 The  onset  of  upwelling  and  development 

of  the  phytoplankton  plume 

It  is  difficult  to  graphically  represent  the  time-de- 
pendent nature  of  the  velocity  field  and  the  corresponding 
spatial  features  of  the  biological  and  chemical  dependent 
variables.  In  the  following  "snapshot"  displays  of  these 
fields,  the  reader  .should  remember  many  functions  (such  as 
growth  of  phytoplankton)  vary  with  time  of  day.  For  example, 
a diel  periodicity  occurs  in  the  concentration  of  ammonia 
and  nitrate  which  is  not  apparent  in  the  following  pictures. 
As  a convention,  the  velocity  field  and  the  corresponding 
spatial  distribution  of  the  dependent  variables  are 


63 


64 


displayed  at  the  end  of  a model  day.  The  nondimens ional 
model  solutions  have  been  multiplied  by  Nt  to  regain  units 
of  concentration  and  thus  help  the  reader  compare  the  model 
results  with  observations. 

The  velocity  field  after  4 days  elapsed  time  in  the 
strong  upwclling  case  is  presented  in  Fig.  18a.  The  linear- 
ized bottom  topography  and  the  wind  stress  at  the  time  of 
the  velocity  field  snapshot  are  also  shown.  The  wind  stress 
vector  has  a magnitude  of  -0.5  dyne  cm  and  points  south- 
ward. The  velocity  field  is  visualized  by  vectors  repre- 
senting the  position  and  instantaneous  velocity  of  tracer 
particles  which  have  been  advectod  by  the  flow.  Only  the 
vectors  within  25  km  of  the  coast  are  shown,  although  the 
circulation  is  calculated  out  to  3100  km.  The  vector  arrows 
are  scaled  by  the  maximum  vector  occurring  in  the  field  at 
that  time.  It  should  be  remembered  that  each  vector's  hori- 
zontal and  vertical  scales  differ  by  two  orders  of  magnitude. 
Two  cyclonically  rotating  circulation  cells  are  evident  in 
Fig.  18a.  The  lower  cell  advects  water  towards  the  coast 
along  the  bottom  and  up  into  the  euphotic  zone  within  5 km 
of  the  coast.  This  upwelled  water  either  continues  to  rise 
to  the  surface  or  moves  offshore  beneath  the  second  counter- 
clockwise rotating  circulation  cell.  In  this  second  cell  the 
flow  is  weakly  onshore  at  30  m depth  and  strongly  offshore 
near  the  surface.  This  two-cell  circulation  is  similar  to 
the  conceptual  diagram  of  coast&l  upwelling  off  Oregon 
presented  in  Fig.  16  of  Mooers , Collins,  and  Smith  (1976). 
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Figure  18b  displays  the  daily  gross  primary  production 
of  the  water  column  calculated  at  each  1 km  grid  point  in 
the  model  region  at  day  4.  The  distribution  of  phytoplank- 
ton within  SO  km  of  the  coast  and  within  50  m of  the  sur- 
face after  4 days  elapsed  time  is  shown  in  Fig.  18c. 

Figure  18  shows  that  a high  *n  phytoplankton  biomass 
develops  near  the  coast  with  the  onset  of  upwclling.  The 
formation  of  a,  plume  is  apparent  in  the  phytoplankton  field. 
The  upwclling  of  water  high  in  nitrate  concentration  into 
the  cuphotic  zone  leads  to  the  increase  in  plant  biomass. 
Even  though  the  growth  rate  of  phytoplankton  is  slower  near 
the.  coast  due  to  the  colder  temperatures  there  (Fig.  16a), 
the  nutrient-limited  plants  bloom  with  the  supply  of  nitrate 
(Fig.  18c). 

After  7 days  the  two- cell  circulation  (Fig.  19a)  and 
the  phytoplankton  plume  (Fig.  19c)  are  well  developed.  A 
maximum  plant  biomass  occurs  at  17  m depth  between  6 and  11 
km  offshore  (Fig.  19c).  This  maximum  has  been  advected  away 
from  the  coast  and  closer  to  the  surface.  The  maximum  also 
occurs  at  a shallower  depth  than  the  steady  state  profile 
maximum  (Fig.  11)  because  of  self-shading  by  the  rapidly 
growing  phytoplankton  (Fig.  19b). 

As  upwelling  at  the  coast  intensifies,  downwelling  de- 
velops in  the  region  between  10  and  25  km  offshore  (Fig. 

4 

20a) . The  contours  of  the  phytoplankton  and  nitrate  fields 
near  the  water  column  bottom  boundary  clearly  show  the 
downwelling  (Figs.  20c  and  21b).  The  phytoplankton  maximum 
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Fig.  21.  (a)  Zooplankton  fiold  after  10  days 
in  the  strong  upvelling  cae'e'.  Contour  intervals  are 
1.8  pgat’  N A-1,  (h)  Nitrate  field  after  10  days  in 
strong, upvelling  case.  Contour  intervals  are  2 vfcat 
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increases  as  it  is  advected  further  offshore  and  closer  to 
the  surface  (Fig.  20c). 

The  zooplankton  field  shows  the  upwelling  and  offshore 
transport  of  water  low  in  herbivore  biomass  (Fig.  21a).  The 
highest  zooplankton  concentrations  occur  seaward  of  the 
phytoplankton  plume. 

The  concentration  of  nitrate  in  the  ouphotic  zone  is 
kept  low  by  plant  production.  Only  where  strong  upwelling 
occurs  does  the  supply  of  NO,  exceed  its  utilization,  allow- 
ing high  concentrations  of  nitrate  to  reach  the  surface 
(Fig.  21b). 

In  the  steady  state  ammonia  profile,  the  maximum  con- 
centration occurs  at  35  m (Fig.  13).  During  the  strong 
upwelling  case,  this  maximum  intensifies  as  nutrient  supply 
to  the  euphotic  zone  leads  to  increased  phytoplankton  pro- 
duction, zooplankton  grazing  and  fecal  pellet  production. 

The  slow  oxidation  of  ammonia  into  nitrite-nitrate  is  ex- 
ceeded by  the  nutrient  input  from  decomposing  fecal  pellets, 
resulting  in  high  concentrations  of  NH^  at  35  to  45  m depth 
(Fig.  22a).  Upwelling  of  water  low  in  ammonia  concentration 
within  10  km  of  the  coast  and  between  33  and  50  km  offshore 
further  intensifies  the  ammonia  gradient. 

Figure  22b  shows  the  development  of  a plume  in  the 

detritus  field.  The  detritus  maximum  at  day  10  in  the  strong 

upwelling  case  occurs  several  meters  below  and  several  kilo- 

% 

meters  seaward  of  the  phytoplankton  plume.  The  main  source 
of  this  detritus  is  fecal  pellets  produced  by  herbivores 


***.  cc.  Kai  Ammonia  field  after  10  days  in  t 
str°ngiupvelling  case.  Contour  intervals  are  0.6  pga 
no11-  SpatiQl  distribution  of  detritus  after 

Rtf0ne  u?weUinS  ^e. . Contour  inter- 
vals are  0.12  ygat*  N £-1. 
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grazing  in  the  phytoplankton  plume.  The  fecal  pellets  sink 
as  they  are  advccted  offshore. 

4.2  Cessation  of  upwelling  and  decay 
of  the  phytoplankton  plume 

When  the  wind  linearly  decrease  to  zero  during  day  11, 
the  upwelling  at  the  coast  begins  to  decline  and  the  offshore 
downwclling  slowly  relaxes  (Figs.  23a,  24a).  The  upwcllcd 
water  is  not  as  strongly  advccted  offshore  and  below  the 
euphotic  zone  as  before.  This  allows  the  phytoplankton  to 
bloom  on  the  nitrate  which  remains  longer  in  the  euphotic 
zone.  The  highest  plant  concentration,  23  pgat  N jf*,  occurs 
on  day  13,  two  days  after  the  wind  stress  falls  to  zero.  As 
upwelling  subsides  the  phytoplankton  maximum  recedes  toward 
the  coast  and  deepens  (Fig.  23c). 

After  20  days  the  plant  maximum  has  decreased  to  IS. 4 
Pgat  N A"1  and  occurs  contiguous  with  the  coast  at  a depth 
of  IS  m (Fig.  24c).  With  no  wind  forcing,  upwelling  ceases 
(Fig.  24a)  and  the  euphotic  zone  becomes  depleted  of  nutri- 
ents (Fig.  25b) . The  zooplankton  standing  stock  declines 
nearshore  (Fig.  25a).  The  detritus  maximum  decays  as  fecal 
pellet  production  lessens  (Fig.  26b).  The  ammonia  maximum 
offshore  begins  to  decrease  at  day  20  (Fig.  26a)  as  more 
NH^  is  oxidized  to  N0j  than  is  added  by  decomposing  fecal 
pellets . 

4 . 3 Daily  primary  production  of  the  water  column 
during  the  strong  upwelUng  case 

The  daily  gross  primary  production  of  the  water  column 
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is  plotted  above  the  corresponding  phytoplankton  field  for 
days  4,  7,  10,  15,  and  20  (Figs.  18b,  19b,  20b,  23b,  and 
24b).  These  figures  show  that  as  upwelling  supplies  nitrate 
to  the  nutrient-limited  cuphotic  zone,  primary  productivity 
increases.  The  plants  grow  fastest  where  the  supply  of 
nutrients  is  greatest.  The  maximum  daily  primary  production 
occurs  on  day  10  (Fig.  20b)  when  upwelling  is  strongest 
(Fig.  20a).  Primary  production  increases  not  only  at  the 
locus  of  upwelling,  but  also  in  surface  waters  which  are 
advcctcd  seaward.  While  the  highest  daily  primary  produc- 
tion occurs  in  the  region  of  most  intense  upwelling  (Figs. 
20a  and  b) , the  maximum  phytoplankton  standing  crop  is  found 
offshore  (Fig.  20c). 

Figure  27  depicts  the  relationship  betwcon  the  north- 
south  component  of  the  wind  stress  and  daily  primary  produc- 
tion. The  increasing  amplitude  of  the  production  curve  in 
Fig.  27  reflects  the  growing  concentration  of  phytoplankton 
biomass  in  the  plume.  Productivity  continues  to  increase 
during  upwelling  conditions,  but  declines  when  the  northerly 
wind  stress  decays  on  day  11.  With  no  further  wind  forcing 
of  upwelling,  phytoplankton  growth  becomes  nutrient- limited. 

Residual  upwelling  accounts  for  the  shoreward  increase 
in  the  productivity  curves  in  Figs.  23b  and  24b.  Residual 
downwelling  in  the  region  10  to  25  km  offshore  reduces  the 
productivity  by  as  much  as  9 mg  N m day  as  plants  are 
advected  to  aphotic  zone  depths.  The  reduction  in  produc- 
tivity in  this  region  is  also  due  to  lower  temperatures. 
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4 . 4 Plume  structure  and  primary  productivity 

during  intermittent  upwelling 

To  investigate  the  response  of  primary  production  to 
intermittent  upwelling,  Thompson's  (1974)  numerical  circula- 
tion model  was  run  with  a variable  wind  stress  calculated 
from  the  August  1973  recordings  of  an  anemometer  located  on 
the  south  jetty  off  Newport,  Oregon.  The  wind  stresses  in 
the  north-south  and  oast-west  directions  are  shown  in  Figs. 
28  and  29,  respectively.  The  circulation  model  was  driven 
from  rest  on  August  1,  1973  with  initial  conditions  speci- 
fied from  observational  hydrographic  data.  The  resulting 
time  dependent  u,  v,  and  w velocities  were  used  to  advect 
the  biological  and  chemical  dependent  variables.  Initial 
conditions  for  the  biological  simulation  model  were  the  same 
as  for  the  strong  upwelling  case,  i.e.,  the  steady  state 
solution  of  the  (x,z,t)  biological-chemical  model  in  the 
absence  of  advection. 

The  two-cell  circulation  found  in  the  strong  upwelling 
case  develops  during  the  first  five  days  of  August  (not 
shown).  The  north-south  wind  stress  during  this  period 
(Fig.  28)  is  similar  in  magnitude  to  the  northerly  wind 
stress  of  the  strong  upwelling  case  (Fig.  27).  A phyto- 
plankton plume  similar  to  that  of  the  strong  upwelling  case 
develops  by  August  5 (Fig.  30a).  Daily  primary  production 
(Fig.  28)  is  comparable  to  that  calculated  for  the  fifth 
day  of  the  strong  upwelling  case  (Fig.  27). 
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Fig.  29.  The  east-vest  component  of  the  wind  stress  calculated 
from  the  Newport,  Oregon,  anemometer  data  recorded  during  August  1973. 
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Fig.  30.  (a)  The  spatial  distribution  of 

phytoplankton  on  August  5 in  the  intermittent  up- 
yelling  case.  Contour  intervals  are  1.2  pgat  N A"1, 
(b)  Same  as  (a)  but  on  August  15.  Contour  intervals 
are  1.4  ygat  N IT1. 
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The  two-cell  circulation  decays  into  a more  confused 
upwelling  circulation  after  August  6 when  periods  of  pro- 
longed relaxation  of  the  wind  stress  occur  (Pig.  28). 

While  primary  production  steadily  increases,  with  continu- 
ous wind  forcing  in  the  strong  upwelling  case  (Pig.  27), 
calculations  of  daily  primary  production  from  August  5 to 
August  20  show  a variability  which  is  Telatcd  to  fluctua- 
tions in  the  north- south  component  of  the  wind  stress 
(Pig.  28).  Each  major  intensification  of  the  longshore 
wind  component  is  followed  by  an  increase  in  primary  pro- 
duction (Figs.  31  and  32).  Upon  relaxation  of  the  wind 
stress,  the  rate  of  upwelling  slows  and  the  primary  pro- 
ductivity decreases.  The  wind  event  of  August  13  to  15 
(Fig.  28)  results  in  advection  of  the  phytoplankton  maximum 
away  from  the  coast  (Pig.  30b),  much  as  observed  after  10 
days  in  the  strong  upwelling  case. 

The  spatial  structure  of  the  dependent  variables  (P, 

Z,  N03,  NH^  and  D)  all  develop  features  during  intermittent 
upwelling  similar  to  that  found  in  the  strong  upwelling 
case.  The  zooplankton,  nitrate,  ammonia  and  detritus  fields 
as  they  appear  on  August  15  are  shown  in  Figs.  33  and  34. 

The  quiescent  period  from  August  19  to  23  is  followed 
by  a reversal  in  the  north-south  wind  component.  Down- 
welling  develops  at  the  coast.  The  contours  of  the  phyto- 
plankton field  on  August  23  (Fig.  35)  clearly  show  the 
downwelling  within  4 km  of  the  coast. 
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Fig*  3^.  The  ammonia  field  on  August.  3$  in  the 
intermittent  upvelling  case.  Contour  intervals  are  0.6 
ygat  HHi t-1.  (b)  The  spatial  distribution  of  detritus 

on  August  15  in  the  interr' ttent  upvelli  g case. 

Contour  intervals  are  0.1  pgat  N l~  . 
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Fig.  35*  The  phytoplankton  field  predicted  by  the  (x,z,t) 
productivity  model  for  August  23  in  the  intermittent  upvelling 
case.  Contour  intervals  are  1.8  ygat  N 
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4 . 5 Model  sensitivity  analysis 

The  (x,z,t)  patchiness  simulation  model  contains  34 
independent  or  specified  parameters  (Table  1).  The  values 
of  many  of  these  parameters  are  chosen  from  observational 
and  laboratory  data  and  the  remaining  are  deduced  from  bio- 
logical and  physical  oceanographic  theory.  It  is  critical 
to  understand  the  response  of  the  model  to  variation  in  the 
parameter  values,  especially  for  those  parameters  whose 
value  is  controversial. 

The  ideal  method  of  parameter  value  investigation  is 
an  analytical  sensitivity  analysis.  However,  the  theory  of 
sensitivity  analysis  is  well  developed  only  for  steady 
state  one-box  models  (Tomovic,  1965).  There  is  an  obvious 
need  to  develop  analytical  tools  to  investigate  the  time 
dependent  nature  of  spatial  models.  For  the  present,  one 
must  resort  to  an  empirical  sensitivity  analysis  for  com- 
plex spatial  models,  such  as  the  one  described  here. 

One  begins  an  empirical  sensitivity  analysis  with  the 
best  estimate  for  all  parameter  values,  and  then  conducts 
a qualitative  investigation  of  the  model’s  response  to 
variation  of  an  individual  parameter.  For  the  sake  of 
brevity  the  results  of  an  extensive  empirical  sensitivity 
analysis  will  be  only  summarized  in  the  following. 

The  spatial  distributions  of  the  dependent  variables 
P,  Z,  N03,  NH^  and  D are  quite  sensitive  to  the  value  of 
the  vertical  diffusion  coefficient,  Ky.  Increased  diffusion 
in  the  vertical  smooths  the  gradient  structure.  With  an 
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increased  value  of  Ky,  the  phytoplankton  plumo  deepens  and 
its  vertical  profile  is  rcorc  smooth.  Nutrients  arc  supplied 
to  the  ouphotic  zone  by  vertical  diffusion  at  a faster  rate 
and  daily  gross  primary  production  increases.  The  spatial 
model  is  loss  responsive  to  an  increase  in  the  horizontal 
diffusion  coefficient,  K^,  although  smoothing  in  the  hori- 
zontal structure  of  all  variables  does  occur. 

The  sinking  rate  of  dead  organic  matter  has  been  re- 
viewed by  Smayda  (1969)  and  more  recently  measured  in  the 
laboratory  by  Smayda  (1970)  and  by  Fowler  and  Small  (1972). 
Depending  on  the  shape  and  density  of  the  particle,  sinking 
rates  range  between  one  and  several  hundred  meters  per  day 
for  dead  phytoplankton  and  fecal  pellets.  Envisioning 
detritus  in  this  model  as  lysed  phytoplankton  cells  and 
floccular, newly  egested  copepod  fecal  pellets,  the  author 
has  used  a detritus  sinking  rate,  w , of  0.01  cm  sec"1  or 
8.6  m day"1  in  all  spatial  solutions  shown  in  Section  4. 

An  increase  in  wg  results  in  a deepening  of  the  maximum  in 
the  steady  state  detritus  profile  and  a more  rapid  loss  of 
limiting  nutrient  from  the  euphotic  zone  in  the  form  of 
decomposing  detritus. 

The  response  of  primary  production  to  variation  in  the 
temperature  field  was  explored  by  assuming  no  offshore  varia 
tion  in  temperature.  With  20°C  surface  temperatures  at  the 
coast,  daily  gross  primary  production  increased  by  65$ 
(compare  Figs.  20b  and  36a).  The  standing  crop  of  phyto- 
plankton at  day  10  in  the  strong  upwelling  case  increased 
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by  2 ygat  N over  the  standing  stock  observod  in  Fig.  20c 
The  phytoplankton  plume  maximum  occurred  closer  to  the  coast 
and  to  the  surface  (Fig.  36b).  Thus  a proper  specification 
of  the  temperature  field  is  necessary  to  prevent  ovcrcstima- 
tion  of  the  primary  production.  Thompson's  (1974)  model 
predicts  the  offshore  variation  in  temperature  of  two  dis- 
tinct, density  layers  off  Oregon.  Given  the  sensitivity  of 
primary  productivity  to  the  temperature  field,  future  model- 
ing efforts  will  incorporate  this  time-dependent  temperature 
prediction. 

Phytoplankton  growth  is  quite  sensitive  to  the  avail- 
ability of  light.  Greater  incident  radiation  at  the  surface 
results  not  only  in  an  enhancement  of  the  light  inhibition 
effect,  but  paradoxically  in  an  increase  in  daily  primary 
production.  The  depth  of  the  cuphotic  zone  dues  not  in- 
crease when  the  dense  phytoplankton  bloom  becomes  self- 
shading. Only  a reduction  in  the  self-shading  coefficient, 
Kp,  or  the  light  extinction  coefficient  for  seawater,  *w, 
would  allow  an  increased  euphotic  zone  depth. 

The  sensitivity  of  the  model  to  the  biological  rate 
processes,  such  as  the  grazing  coefficient,  Rm,  is  best  de- 
termined by  an  analytical  sensitivity  analysis  of  the 
one-box,  nonspatial  nondimensional  model.  Empirical  anal- 
ysis of  (13)  - (17)  show  the  parameters  describing  the 
zooplankton  dynamics;  i.e.,  grazing,  egestion  and  excretion 
are  most  important  relative  to  the  phytoplankton  growth  rate, 
V , in  determining  the  steady  state  concentration  of  the 
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dependent  variables  P,  Z,  NO-,  NH^  and  D.  O'Brien  and 
Wroblcwski  (1976)  showed  conclusively  phytoplankton  growth 
and  nutrient  recycling  by  zooplankton  are  the  most  impor- 
tant and  crucial  processes  in  their  (x,z,t)  spatial  eco- 
system model  for  the  west  Florida  continental  shelf. 

Almost  the  same  steady  state  solution  of  nonspatial 
(13)  - (17)  is  attained  using  a zooplankton  egestion  formu- 
lation as  recommended  by  Conover  (1966)  where  egestion  is 
a constant  fraction  of  ingestion.  In  the  (z,t)  or  (x,z,t) 
model  where  small  concentrations  of  phytoplankton  occur, 
the  choice  of  formulation  does  lead  to  important  differences 
in  egestion  rate. 

The  detritus  regeneration  rate  parameter,  «t> , and  the 
ammonia  oxidation  rate  parameter,  ft,  have  long  o-folding 
time  scales,  yet  their  values  arc  important  in  the  time- 
dependent,  spatial  solutions.  Increasing  the  rate  of 
decomposition  of  detritus  into  ammonia,  for  example,  libe- 
rates more  nutrient  for  plant  growth  in  the  euphotic  zone 
and  results  in  a higher  plant  standing  crop.  Increasing 
the  rate  of  oxidation  of  ammonia  can  prevent  the  development 
of  an  ammonia  maximum  in  the  spatial  solutions.  The  impor- 
tance of  the  small  valued,  senescent  cell  lysis  coefficient, 
S,  becomes  evident  in  the  (z,t)  and  (x,z,£)  solutions,  where 
this  term  leads  to  the  low  concentrations  of  pli-nts  in  the 
aphotic  zone. 

A positive  grazing  threshold,  P^,  is  not  necessary 
to  prevent  phytoplankton  from  being  locally  grazed  to 
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extinction.  Diffusion  of  plants  from  areas  of  high  concen- 
tration to  low  concentration  keeps  a finite  amount  of  phyto- 
plankton standing  crop  present  at  each  grid  point.  This  is 
an  important  result.  In  essence^  is  a parameterization  of 
turbulent  diffusion  in  nonspatial  models  (Wroblewski  and 
O'Brien,  1976). 

4 . 6 Comparison  with  observations 

To  establish  credibility  of  the  model,  the  spatial 
solutions  in  Sections  4.1  to  4.4  are  compared  with  actual 
observations.  It  should  be  remembered  that  this  model  is 
an  idealized  description  of  the  complex  upwelling  ecosys- 
tem whose  irregular  physical  boundaries  are  three-dimen- 
sional. Nevertheless,  if  tho  model  is  to  be  a pedagogic 
tool,  its  major  features  must  be  verifiable  by  observations. 

Many  of  the  features  found  in  the  model  solutions  have 
been  observed  off  the  Oregon  coast  during  the  Coastal  Up- 
welling  Experiment  or  have  been  described  in  the  literature. 
Fig.  37  shows  the  development  of  a phytoplankton  plume  off 
Newport,  Oregon,  in  August  1972.  The  similarity  in  orienta- 
tion, position  and  length  scale  between  the  simulated  and 
observed  plume  structure  is  striking. 

Turbidity  measurements  made  by  Pak  et  al.  (1970)  in 
June  1967,  indicate  a suspended  particle  maximum  existed 
within  the  same  depth  range  and  offshore  location  as  the 
phytoplankton  and  detritus  maxima  shown  in  Figs.  30b  and 
34b.  To  explain  this  particle  distribution,  they  proposed 
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a sequence  of  events  in  which  nutrient  rich  water  upwells 
close  to  shore,  becomes  increasingly  concentrated  with 
particles  of  biological  origin  as  it  moves  offshore,  and 
gradually  sinks  below  the  permanent  pycnoclinc  which  slopes 
up  to  near  the  surface.  This  mechanism  was  adopted  by 
Small  and  Ramberg  (1971)  and  Anderson  (1972)  to  explain  the 
plume  feature  often  observed  in  the  phytoplankton  distribu- 
tion during  upwelling.  Kitchen  ct  al . , (1975)  have  shown 
the  particle  maximum  in  July  1973,  was  composed  mostly  of 
particles  20  to  50  ym  in  diameter  nearshore  and  of  increas- 
ingly smaller  sized  particles  with  progression  offshore. 

They  concluded  the  suspended  particles  in  the  upwelling 
region  close  to  the  coast  were  diatoms,  which  were  progress 
sively  replaced  by  dinoflagellates  in  the  offshore  nutrient 
limited  waters. 

Peterson  (1972)  observed  the  zooplankton  standing 
stock  was  greater  in  the  oceanic  region  than  on  the  conti- 
nental shelf  off  Oregon  in  the  summer.  Only  during  winter 
was  there  a high  zooplankton  standing  stock  close  inshore, 
where  the  primary  production  is  greatest.  Peterson  sug- 
gested predation  upon  zooplankton  nauplii  and  offshore 
water  transport  as  reasons  for  low  coastal  herbivore  bio- 
mass. The  model  solutions  indicate  offshore  transport 
leads  to  a reduction  of  zooplankton  standing  stock  over  the 
shelf  (Figs.  25a  and  33a).  Advective  transport  could 
explain  the  observed  high  coastal  zooplankton  standing 
stock  in  winter  when  the  Ekman  surface  flow  is  Dredominantly 


I 


08 

onshore  (Bakun,  1973). 

Peterson  and  Miller  (1975)  found  high  zooplankton 
biomasses  occurred  within  10  km  of  the  coast  during  the 
1969-71  upwolling  seasons,  as  did  Myers  (1975)  in  August 
1S73.  Standing  stocks  were  highest  within  2 km  of  the  beach, 
indicating  the  source  of  these  animals  was  the  coastal  popu- 
lations. If  the  (x,z,t)  model  had  included  in  its  coastal 
boundary  conditions  the  inshore  source  of  zooplankton,  or 
if  zooplankton  were  allowed  to  vertically  migrate  to  depths 
where  onshore  flow  occurs,  the  observations  of  Peterson, 
Miller,  and  Myers  might  have  been  more  closely  reproduced. 

Measurements  of  the  spatial  distribution  of  nitrate  by 
Ball  (1970),  Atlas  (1973),  and  Myers  (1975)  indicate  upwcll- 
ing  of  water  high  in  nitrate  concentration  occurs  at  the 
coast.  Comparison  of  these  observed  nutrient  fields  with 
Figs.  21b  and  33b  shows  good  qualitative  agreement.  Few 
measurements  of  the  ammonia  field  were  made  before  the 
Coastal  Upwelling  Experiment.  The  ammonia  data  from  1972 
and  1973  docs  indicate  a maximum  between  10  and  20  m depth 
seaward  of  10  km.  Not  enough  profiles  were  made,  however, 
to  validate  the  offshore  ammonia  maximum  predicted  in  the 
model  solutions  (Figs.  22a  and  34a).  The  NH^  maximum  is 
a temporary  feature  which  decays  when  nutrient  input  by 
upwelling  ceases  and  thus  may  easily  be  missed  in  field 
sampling.  The  spatial  structure  of  ammonia  should  be  more 
thoroughly  investigated  in  future  observational  programs. 
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Field  measurements  of  net  primary  productivity  by  the 

Cl«  technique  performed  by  Anderson  (1964)  and  by  Small  ct 

al,  (1972)  off  Oregon  during  upwelling  ranged  from  O.S  to 
- 2 - 1 

1.5  g m day  . If  a C/H  ratio  of  7 is  assumed  (Small 

and  Rambcrg,  1971)  the  observed  range  in  terms  of  nitrogen 

is  71  to  214  mg  N m day  . Gross  primary  productivities 

calculated  for  the  model  solutions  range  between  78  and  296 

mg  N m day  during  the  strong  upwelling  case,  and  from 

-2  -1 

78  to  226  mg  N m day  during  the  intermittent  upwelling 
case.  This  is  excellent  agreement  between  observation  and 
model.  The  predicted  increase  in  daily  gross  primary  pro- 
duction of  the  water  column  after  a strengthening  of  north- 
erly winds  for  several  days  during  the  intermittent  upwcll- 
ing  case  (Fig.  28)  is  suggested  in  the  CUE  II  productivity 
data  (Small,  in  preparation).  The  bloom  of  phytoplankton 
upon  relaxation  of  winds  after  a major  wind  event  (Section 
4.2)  is  a model  feature  clearly  observed  during  the  CUE 
field  study. 


5.  SUMMARY , CONCLUSIONS  AND  CRITIQUE 


Equations  describing  the  mesoscalc  distribution  and 
ecosystem  dynamics  of  phytoplankton,  herbivores,  detritus 
and  the  nutrients,  nitrate  and  ammonia,  have  been  formu- 
lated for  a time-dependent,  two-dimensional,  numerical 
model  of  phytoplankton  patchiness  during  Oregon  coastal 
upwclling.  New  formulations  for  herbivore  egestion  and 
ammonia- inhibition  of  nitrate  uptake  by  phytoplankton  are 
presented. 

The  goal  of  this  research  has  been  to  provide  a dy- 
namical explanation  for  the  spatial  features  which  are 
consistently  observed  in  the  phytoplankton,  nitrate  and 
detritus  fields  during  upwclling  conditions  off  Oregon. 
While  it  is  not  implied  that  conclusive  verification  of 
the  model  lies  in  the  cited  observations,  the  similarity 
between  model  solutions  and  field  data  suggests  the  complex 
(x,z,t)  model  may  include  the  basic  biological,  chemical 
and  physical  dynamics  governing  primary  productivity  in  the 
Oregon  upwelling  ecosystem. 

Production  is  a function  of  the  availavility  of  light 
and  nutrients,  and  its  rate  is  governed  by  temperature. 
Because  the  physical  and  chemical  environment  in  a coastal 
upwelling  region  is  highly  variable,  production  must,  be 
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modeled  within  a temporal,  spatial  framework.  Smayda  (1966) 
had  only  limited  success  in  predicting  daily  primary  produc- 
tion from  an  empirical  equation  relating  phytoplankton  bio- 
mass directly  to  the  wind  stress  and  surface  tomperaturo. 
Small  et  aU , (1972)  concluded  there  was  no  simple  empiri- 
cal relationship  between  daily  primary  production  off 
Oregon  and  such  environmental  factors  as  incident  radiation, 
photic  depth  or  rate  of  upwclling.  Iverson  et^  al.  , (1974) 
were  able  to  model  phytoplankton  standing  crops  observed  in 
Auke  Bay,  Alaska,  after  allowing  for  the  time-dependent 
supply  of  nitrate  to  the  cuphotic  zone  by  wind-mixing. 

Walsh  (1975)  simulated  phytoplankton  and  nitrate  concentra- 
tions observed  in  the  Peruvian  upwelling  region  with  a 
spatial,  although  time-invariant,  ecosystem  model. 

A detailed  mechanistic  model  is  required  to  give  ade- 
quate predictions  of  primary  production.  Nonspatial  eco- 
system models  such  as  that  of  Cushing  (1971)  are  useful 
for  investigating  the  relationship  between  productivity 
and  lighfej  grazing  and  rate  of  upwelling,  but  they  neglect 
the  important  influence  of  a variable  environment  upon 
primary  productivity. 

The  dominant  role  cf  advoction  in  determining  the 
spatial  configuration  of  plankton  and  nutrient  fields  in 
upwelling  regions  is  obvious  from  the  model  solutions.  To 
properly  delineate  the  mesoscale  phytoplankton  plume  off 
Oregon,  a 2.5  m vertical  and  1 km  horizontal  resolution  of 
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the  velocity  field  was  required.  Finer  grids  arc  necessary 
to  model  smaller  scale  phytoplankton  patchiness. 

The  physical  mechanisms  which  lead  to  the  phytoplank- 
ton plume  structure  observed  over  the  continental  shelf 
off  Oregon  appear  to  be  those  suggested  by  Pak  et  al. , (1970) 
and  described  dynamically  by  Mooers  e£  al . , (1976).  The 
numerical  upwclling  model  of  Thompson  (1974)  simulates  this 
circulation  i.c.;  upwclling  at  the  coast,  offshore  transport 
and  downVclling  in  the  region  10  to  25  km  offshore.  This 
physical  model  when  coupled  to  the  biological  and  chemical 

dynamics  described  in  Section  2,  reproduces  the  phytoplankton 

\ 

plume  often  observed  off  Oregon  during  the  upwelling  season. 
However,  the  two-cell  circulation  is  but  one  of  several  up- 
wclling patterns  believed  created  by  strong  northerly  winds 
off  Oregon. 

In  spite  of  the  successful  reproduction  of  the  biologi- 
cal and  chemical  features  observed  off  Oregon  during  upwelling, 
this  model  is  fundamentally  limited  in  its  ability  to  predict 
the  distribution  of  its  dependent  variables  in  tine  and  spare. 
Numerical  model  solutions  of  a nonlinear,  dissipative  open 
thermodynamic  system  (such  as  the  Oregon  upwelling  system) 
will  diverge  from  the  "real"  solution  if  the  initial  state 
of  the  system  is  not  specified  exactly  (Lorenz,  1969).  In 
other  words,  one  must  know  the  actual  initial  conditions  for 
all  dependent  variables  for  all  time  and  space.  The  inevit- 
able discrepancy  between  the  model  prediction  and  observa- 
tions stems  from  the  fact  that  continuous  physical  processes 
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in  the  ocean  are  represented  by  finite  approximations. 

Platt  et  al.  , (1975)  estimate  the  limit  to  predictability 
of  features  with  length  scales  of  the  order  10  km  in  the 
ocean  is  at  best  several  days.  After  this  period,  the 
spatial  solution  and  the  observed  fields  may  differ  appreci- 
ably. 

An  important  criticism  of  this  model  is  the  assumption 
that  all  parameter  values  have  no  spatial  or  temporal  vari- 
ability. For  example,  the  light  extinction  coefficient  for 
the  seawater  off  Oregon  was  assumed  constant  with  distance 
offshore.  The  extinction  coefficient  is  known  to  decrease 
with  distance  from  shore  as  turbidity  caused  by  particles 
suspended  in  river  runoff  decreases  (Small  and  Curl,  1968). 
If  a spatially  variable  k had  been  assumed,  the  subsurface 

W 

phytoplankton  maximum  observed  by  Anderson  (1969)  in  oceanic 
waters  off  Oregon  may  have  been  reproduced  close  to  the  sea- 
ward boundary  cf  the  model.  Instead,  this  feature  is  missed 
altogether.  Nevertheless,  the  model's  spatial  solutions  in 
the  region  of  interest  arc  reasonable. 

Longshore  variation  in  all  model  variables  has  been 
ignored,  although  Peffley  and  O'Brien  (1976)  have  demon- 
strated conservation  of  mass  in  the  (x,z)  plane  is  not 
achieved  during  Oregon  upwelling.  However,  high  altitude 
photographs  (Keen  and  Pearcy,  1973)  and  airborne  remote 
sensing  of  sea  surface  temperature  (Holladay  and  O'Brien, 
1975)  show  sea  surface  colors  and  temperature  isotherms 
tend  to  be  parallel  to  the  coast  for  several  tens  of 
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kilometers.  This  suggests  there  is  less  longshore  variabil- 
ity in  the  intensity  of  upwclling  off  the  Oregon  coast  than 
observed  in  other  upwclling  regions  - notably  Peru,  where 
phytoplankton  are  distributed  in  shallow  tongues  (Kelley  et 
al. , 197S) . The  model  also  ignores  the  possible  influence 
of  the  Columbia  River  effluent  tongue  upon  the  physical  and 
chemical  environment  off  Oregon. 

This  model  represents  the  most  complex,  time-dependent, 
spatial,  numerical  model  of  an  upwelling  ecosystem  developed. 
A major  strength  of  the  model  is  the  formulation  of  numerous 
physical,  chemical  and  biological  processes  into  an  inte- 
grated framework.  A major  weakness  is  tho  treatment  of  in- 
dividual phytoplankton  and  zooplankton  species  as  members 
of  a food  web  with  invariable  tropho- dynamic  relationships. 
Few  organisms  play  such  a simple  structured  role  in  nature. 


APPENDIX 


Definition  of  Symbols  and 
Scaling  Relationships 


Dimensional 

Quantity 

Definition 

Scaling 

Factor 

Nondimensional 

Quantity 

- 

Temperature  function 
parameters 

•- 

a,b 

c 

Temperature  function 
parameter 

- 

- 

D 

Detrital  nitrogen 

D/Nt 

D' 

d 

Daylcngth  fraction  of 
a day 

- 

- 

E 

m 

Maximum  herbivore 
egestion  coefficient 

VVm 

P 

H 

Characteristic  verti- 
cal length  scale 

- 

- 

rs 

Light  saturation  para- 
meter when  8 * 0 or 
n * 0 

T 

- 

Light  intensity  at 
local  apparent  noon 

- 

- 

Kh 

Horizontal  eddy 
diffusivity 

- 

- 

Kv 

Vertical  eddy 
diffusivity 

- 

- 

Half -saturation  con- 
stant for  growth 

- 

- 

ku 

Half-saturation 
constant  for  nutrient 
uptake 

VNt 

a 
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Dimensional 
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Definition 

Scaling 

Nondimcnsional 

Quantity 

Factor 

Quantity 

L 

Characteristic  horizontal 

length  scale 

- 

• 

nii4 

Ammonia 

NH4/Nt 

my 

NO- 

0 

Nitrate  plus  nitrite 

N03/Nt 

no3- 

Nt 

Total  concentration  of 
biologically  limiting 
nutrient 

• 

p 

Phytoplankton  nitrogen 

P/Nt 

p' 

pt 

Herbivore  grazing 
threshold 

Vt 

p* 

Herbivore  maximum 
grazing  rate 

b 

T 

Temperature 

- 

- 

t 

Time 

tvm 

T 

u 

x-directed  velocity 
component 

u/U 

U* 

U 

Typical  value  of  the 
horizontal  velocity 

- 

- 

Vm 

Phytoplankton  maximum 
nutrient  uptake  rate 

- 

- 

W 

Typical  value  of  the 
vertical  velocity 

- 

- 

w 

z-directed  velocity 
component 

w/W 

w' 

ws 

Sinking  velocity  of 
detritus 

vs/K 

V 

X 

Tangent-plane  Cartesian 
coordinates:  x positive 
toward  the  coast 

x/L 

x' 

Z 

Zooplankton  nitrogen 

Z/Nt 

Z’ 

z 

Tangent-plane  Cartesian 
coordinates:  z positive 

upward 

z/H 

z’ 

I 

I 

1 

1 

I 

I 

I 

I 

I 

l 

l 

I 

I 

I 

I 

f 

( 


i 


107 


Dimensional 

Quantity 


Definition 


r Herbivore  excretion 
coefficient 


Scaling  Nondimensional 
Factor  Quantity 


r /v 

m 


Y 


A Herbivore  ogestion  rate 

at  phytoplankton  concen- 
tration P«.  A/V„ 

t m 

Light  function  parameters  - 

k Light  extinction  per  unit 

" phytoplankton  nitrogen 

k Light  extinction  cocffi- 

w cient  of  pure  seawater 

and  any  nonphytoplanktonic 
material 

A Ivlev  constant  ANt 

u Phytoplankton  maximum 

m growth  rote 

H Phytoplankton  nutrient 

loss  coefficient  E/Vm 

Pi 

T Herbivore  egestion  T 

coefficient 

<J>  Detritus  decomposition 

parameter  'fr/V^ 

¥ Nitrate  uptake  inhibi- 
tion parameter  ¥Nt 

ft  Ammonia  oxidation 

coefficient  ft/Vm 


6 


n,G 


X 


5 

TT 

U 


<!* 

* 
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